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The phase diagram and critical properties of the A^-component London superconductor are stud- 
ied both analytically and through large-scale Monte-Carlo simulations in d = 2 -I- 1 dimensions 
(components here refer to different replicas of the complex scalar field). Examples are given of 
physical systems to which this model is applicable. The model with different bare phase stiffnesses 
for each component, is a model of superconductivity which should arise out of metallic phases of 
light atoms under extreme pressure. A projected mixture of electronic and protonic condensates in 
liquid metallic hydrogen under extreme pressure is the simplest example, corresponding to A*' = 2. 
These are such that Josephson coupling between different matter field components is precisely zero 
on symmetry grounds. The A''-component London model is dualized to a theory involving vortex 
fields with highly nontrivial interactions. We compute critical exponents a and u for N = 2 and 
A'^ = 3. Direct and dual gauge field correlators for general A'' are given and the N = 2 case is studied 
in detail. The model with N = 2 shows two anomalies in the specific heat when the bare phase 
stiffnesses of each matter field species are different. One anomaly corresponds to an inverted SDxy 
fixed point, while the other corresponds to a 3Dxy fixed point. Correspondingly, for A^ = 3, we 
demonstrate the existence of two neutral SDxy fixed points and one inverted charged SDxy fixed 
point. For the general case, there are A'^ fixed points, namely one charged inverted SDxy fixed point, 
and N — 1 neutral SDxy fixed points. We explicitly identify one charged vortex mode and A'^ — 1 
neutral vortex modes. The model for N — 2 and equal bare phase stiffnesses corresponds to a field 
theoretical description of an easy-plane quantum antiferromagnet. In this case, the critical expo- 
nents are computed and found to be non SDxy values. The A''-component London superconductor 
model in an external magnetic field, with no inter-species Josephson coupling, will be shown to have 
a novel feature, namely A'^ — f superfluid phases arising out of iV charged condensates. In particular, 
for N = 2we point out the possibility of two novel types of field-induced phase transitions in ordered 
quantum fluids: i) A phase transition from a superconductor to a superfluid or vice versa, driven 
by tuning an external magnetic field. This identifies the superconducting phase of liquid metallic 
hydrogen as a novel quantum fluid, ii) A phase transition corresponding to a quantum fluid analogue 
of sublattice melting, where a composite field-induced Abrikosov vortex lattice is decomposed and 
disorders the phases of the constituent condensate with lowest bare phase stiffness. Both transitions 
belong to the SDxy universality class. For N > S, there is a new feature not present in the cases 
A'" = 1 and N = 2, namely a partial decomposition of composite field-induced vortices driven by 
thermal fluctuations. A "color electric charge" concept, useful for establishing the character of these 
phase transitions, is introduced. 
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INTRODUCTION 



Ginzburg-Landau (GL) theories with several complex 
scalar matter fields minimally coupled to one gauge field 
are of interest in a wide variety of condensed matter sys- 
tems and beyond. This includes such apparently dis- 
parate systems as the two-Higgs doublet modeli, super- 
conducting low temperature phases of light atoms such 
as hydrogenSi^ under extreme enough pressures to pro- 
duce liquid metallic states, and effective theories for easy- 
plane quantum antiferromagnets-'-^'^. Well known cases 
of multicomponent systems are represented by multiband 
superconductors'" like MgB2 where there are two or- 
der parameters corresponding to Cooper pairs made up 
of electrons living on different sheets of Fermi surface. 
In that case however condensates are not independently 
conserved and the U{1) x t/(l) symmetry is broken to 



[/(I), so the main results of this paper do not apply to 
multiband superconductors. In contrast, in the projected 
liquid metallic state of hydrogenSi^, which appears being 
close to a realization in high pressure experimentsAS, the 
scalar fields represent Cooper pairs of electrons and pro- 
tons. This excludes, on symmetry grounds, the possibil- 
ity of inter-flavor pair tunneling, i.e. there is no intrin- 
sic Josephson coupling between different species of the 
condensate. This sets it apart from systems with multi- 
flavor electronic condensates arising out of superconduct- 
ing order parameters originating on multiple-sheet Fermi 
surfaces, such as is the case in MgB2. For the latter 
system, Josephson coupling in internal order parameter 
space cannot be ruled out on symmetry grounds, and 
must therefore be included in the description. This is so 
because the Josephson coupling represents a singular per- 
turbation and can never be ignored on sufficiently long 
length scales. This is otherwise well known from stud- 
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ies of extremely layered superconductorsiS, where the 
critical sector is that of the 2'Dxy model in the absence 
of Josephson coupling, while any amount of interlayer 
phase-coupling (in an extended system) produces a crit- 
ical sector belonging to the SDxy universality class. It is 
precisely the lack of Josephson coupling in certain, but by 
no means all systems with multiple flavor order param- 
eters, that opens up the possibility of novel and inter- 
esting critical phenomena. However, even in inter-band 
Josephson-coupled condensates, interesting physics arises 
at finite length scalesii^i^. 

A two-component action with no Josephson coupling 
in (2 -|- 1) dimensions, with matter fields originating in a 
bosonic representation of spin operators, is also claimed 
to be the critical sector of a field theory separating a 
Neel state and a paramagnetic (valence bond ordered) 
state of a two dimensional quantum antiferromagnet at 
zero temperature with easy-plane anisotropy^'^. This 
happens because, although the effective description of 
the antiferromagnet involves an a priori compact gauge 
field, it must be supplemented by Berry-phase terms in 
order to properly describe S = 1/2 spin systema^^ii^. 
Berry-phase terms in turn cancel the effects of monopoles 
at the critical poinli^. Hence, an effective description 
in terms of two complex scalar matter fields coupled to 
one non-compact gauge field suffices to describe the non- 
trivial quantum critical point separating a state with 
broken internal SU{2) symmetry and a paramagnetic 
5L/(2)-symmetric state with broken external symmetry 
(lattice translational invariance). The latter state is the 
valence-bond ordered state. Critical behavior separat- 
ing states differing in this manner is not captured by 
the Landau- Wilson-Ginzburg paradigm^*i^, and requires 
a description of a phase transition without a local or- 
der parameter. An example of such a description is the 
well known Kosterlitz-Thouless phase transition taking 
place in the 2Dxy modeli^. The difference from the 
Kosterlitz-Thouless case and the quantum critical behav- 
ior described above is that while the low temperature 
phase of the 2'Dxy model is a Gaussian fixed line, this 
is not so for either side of the quantum critical point of 
the easy-plane quantum antiferromagnet^i^*Sii^. We also 
mention that another example of a multicomponent sys- 
tem with no inter-component Josephson effect are spin- 
triplet superconductors which are well known to allow 
a variety of topological defects and phase transitionsi2i. 
Some of the topics we discuss below are related to the 
models of spin-triplet paired electrons-^. 

Since the condensates described above in the context 
of light atoms and easy-plane quantum antiferromagnets 
are gauge-charged condensates, the order parameter fla- 
vors are all coupled to each other via a non-compact 
gauge field. This coupling is vastly different from the 
Josephson coupling in the sense that while an iV-flavor 
order parameter condensate with no coupling between 
different species in general will have N phase transi- 
tions, a Josephson coupling between a pair of order pa- 
rameter species will collapse the two independent phase 



transitions they undergo with no coupling, down to one. 
Josephson coupling between all pairs of order parameter 
species will collapse all N phase transitions down to a 
single one, namely an inverted SDxy transition. On the 
other hand, N order parameter species coupled to one 
and the same gauge field will still undergo in general N 
phase transitions, namely one inverted SDxy transition 
where a Higgs phenomenon takes place, followed by N— 1 
SDxy transitions as the coupling constants are increased 
beyond the Higgs/SDxy critical point^iiSi. 

A special feature is presented by the important case 
N ^ 2. Here, it turns out that the dual description of 
the theory is isomorphic to the starting poinlj^iSiiS. Nor- 
mally, in d = 2 -|- 1, a gauge theory dualizes into a global 
theory and vice versa. In contrast a C/(l) x U{1) gauge 
theory dualizes into another U{1) x U (1) gauge theory, i.e. 
the theories are self-dual. In general the theory has two 
separate critical points, one inverted SDxy and one SDxy 
critical point For the special case where the bare phase 
stiffnesses of the two matter fields are equal, as they nat- 
urally are in the case of easy-plane quantum antiferro- 
magnets in the absence of an external magnetic field^, 
another interesting feature appears. In this case, there is 
only one critical point separating two phases described by 
self-dual field theories. This cannot be either an inverted 
SDxy or a SDxy fixed point. Self-duality also precludes 
the possibility of a Z{2) universality class although the 
exponent v that we find for this case appears to be close 
to the Ising value (while a is not). This phase transition 
therefore defines a new universality class, namely that of 
the d = 2 + 1 U{1) X U{1) self-dual gauge theory. 

What happens to such multi-component charged con- 
densates in three dimensions in the absence of Josephson 
coupling between the order parameter components, but 
in the presence of an external magnetic field, has been 
recently studied in Rcfs. 20 21 for the case N = 2, with 
particular emphasis on applications to liquid metallic hy- 
drogen. In this paper, we extend on this and consider in 
detail the effects of tuning the external magnetic field 
and temperature when also TV > 3. New features appear 
compared to the N = 2 case, because composite vortices 
consisting of non-trivial windings in all order parame- 
ter components can now undergo partial decompositions 
by tearing vortices of individual order parameter com- 
ponents off the composite vortices, one after the other. 
We provide a dual picture of these processes: i) as a, vor- 
tex loop proliferation in the background of a composite 
vortex lattice, and ii) as a. metal-insulator transition in a 
system consisting of several "colors of electric charges" in 
a multi color dielectric background. The new concept of 
"color charge" will be introduced and explained in detail 
in this paper. It allows us to determine the universality 
class, and the partially broken symmetries of the partial 
decomposition transitions taking place in multi flavor su- 
perconductors in an external magnetic field. We also 
show that the number of colors A'coior of dual charges ex- 
ceeds the number of field components (flavors) for TV > 3. 

The outline of the paper is as follows. The first six 
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sections of the paper deal with results in zero external 
magnetic field. In Sections VII and VIII we present re- 
sults in finite magnetic field. Readers who wish to consult 
results on finite magnetic field may proceed directly to 
Section VII. 



In Section ^ we introduce the model and the main 
approximation we will use to study the model, as well as 
the duality transform that will be used extensively, along 
with the explicit vortex representation of the model. In 
Section IIIII we explicitly transform the action for the 
N — 2 case into an action consisting of two parts: i) 
one charged vortex mode with vortex interactions medi- 
ated by a massive vector field, and ii) one neutral vortex 
mode with vortex interactions mediated by a gauge field. 
In Section IIVI we compute gauge field correlators and 
dual gauge field correlators in terms of vortex correla- 
tors. This explicitly identifies the mechanism by which 
a thermally driven vortex loop proliferation destroys the 
Higgs phase (Meissner effect) and dual Higgs phasei^^. 
Gauge field correlators are useful in characterizing the 
charged fixed point of the iV-flavor London modeliSiS^, 
while dual gauge field correlators are also useful in char- 
acterizing the N — 1 neutral fixed points^S. In Section 
rVT we present large-scale Monte Carlo (MC) simulations 
for the case N = 2, computing critical exponents at the 
neutral and charged fixed points, as well as the mass of 
the gauge field as a function of temperature. The neu- 
tral fixed point is found to be in the 3Dxy universality 
class, while the charged fixed point is shown to be in 
the inverted 3T)xy universality class. We also consider 
in detail the case when the two bare phase stiffnesses of 
the model are identical, showing that the resulting one 
fixed point is in a new universality class distinct from 
the 3Dxy and inverted 3Dxy universality classes. In Sec- 
tion IVII we present corresponding results for the case 

= 3. In Section IVIll we outline the phases to expect 
for the case N = 2 when an external magnetic field is 
applied. We also present results from large-scale Monte- 
Carlo simulations revealing a novel phase transition in 
the SDxy universality class inside the Abrikosov vortex 
lattice phase at low magnetic fields when temperature is 
increased. 



II. 



MODEL AND DUAL ACTION 



In Section rVIIII we do the same when N > 2, 
emphasizing the qualitatively new features compared to 
the case N = 2. We also introduce a useful "color charge" 
picture of the various partial decomposition transitions 
of the composite vortex lattice that we encounter for the 
case when > 3. In Section ^ we summarize our re- 
sults. In Appendix ^ we identify charged an neutral 
vortex modes for general N. In Appendix ^ we derive 
the vortex representation for the general- A^ case. In Ap- 
pendices El and ^ we derive expressions for gauge field 
correlators and dual gauge field correlators, respectively. 
In Appendix IeI we generalize our dual representation for 
arbitrary N to also include inter-flavor Josephson cou- 
pling. In Appendix ^ we consider Kosterlitz-Thouless 
transitions for the general-A case in two spatial dimen- 
sions at finite temperature. 



For an analysis of the possible phase transitions in 
a GL model of N individually conserved bosonic mat- 
ter fields, each coupled to one and the same t/(l) non- 
compact gauge field, we study a version of the iV-fiavor 
GL theory in (2 + 1) dimensions with no Josephson cou- 
pling terms between order parameter components. More- 
over, we ignore mixed gradient terms, such that there is 
no Andreev-Bashkin effecliS^. The model is defined by N 
complex scalar fields {"i''^\r) \ a = 1,...,A} coupled 
through the charge e to a fluctuating gauge field A(r), 
with the action 



■ N 

E 

.a=l 



|(V-ieA(r))<^(r 



+ nK"'(r)}) + i(VxA(r))^ 



(1) 



where M*^"^ is the mass of the condensate species a. As- 
suming that the individual condensates are conserved, 
the potential must be function of |^'[)"^(r)|2 

only. In this paper, we focus on the critical phenomena 
and phase diagram of Eq. in zero as well as finite ex- 
ternal magnetic field, and for these purposes the model 
in Eq. ^ will be studied in the phase-only approxima- 
tion *[)"^(r) = |4'[,"^|exp[i6'(")(i")] where |*[,"^| is a con- 
stant, i.e. we freeze out amplitude fluctuations of each 
individual matter field. The model we study is therefore 
the generalization to arbitrary N of the frozen-amplitude 
one gap lattice superconductor model also known as the 
London superconductor model^. 

One may well ask what confidence one should put in 
the phase only approximation for all fields when the bare 
phase stiffness of each individual condensate is very dif- 
ferent, such as is the case in LMH. The answer is that one 
can be quite confident that this is a useful and reasonable 
approximation. Consider first the case N = 2. We use 
the phase only approximation with confidence for consid- 
ering the criticality here. It certainly works at the lowest 
critical temperature. After that point, we are left with 
a one-component superconductor. What the field with 
the lowest phase stiffness does above the lowest critical 
temperature is not of interest, it is only the remaining 
field with criticality at higher temperature that matters. 
Hence, significantly above the lowest critical tempera- 
ture, we may still apply the phase only approximation 
for the remaining one-component case. For this field, we 
may use the phase only approximation up to the high- 
est critical temperature with the same confidence as we 
can use the phase only approximation for the field with 
the lowest phase stiffness up to and slightly above the 
lowest critical temperature. The same argument can be 
repeated for arbitrary N: We can use the phase only ap- 
proximation for the fields up to and slightly above their 
respective critical temperatures. After that it is imma- 
terial what they do, it is only the remaining components 
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that matter. 



B. Separation of variables 



A. Basic properties of the model 

Varying Eq. with respect to A, we obtain the equa- 
tion for the supercurrent 

N 



-2e 



2 / 1*0 



(")|2' 



(2) 



Vortex excitations in such an A^-flavor GL model carry 
fractional flux. Consider a vortex where the phase 0'^''^ (r) 
has a 2tt winding around a vortex core, while other phases 
do not have nontrivial windings. Expressing A from Eq. 

, and integrating along a path around the vortex core 
at a distance larger than the magnetic penetration length, 
we obtain an expression for the magnetic flux encom- 
passed by the path given by 



= i Ad\ = 



I* 



(v)\2 
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Q = l 
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(3) 



where $o = 2.07 • lO^^^Tm^ is the flux quantum. As it 
will be clear from a discussion following Eq. (see 
Eq. H12|l 'l. such a vortex has a logarithmically divergent 
energy^'^'^^. Only a composite vortex where all phases 
^("^ have 27rn winding around the core carries integer flux 
and has finite energy. As detailed below, the composite 
vortices are responsible for the magnetic properties of the 
system at low temperatures while thermal excitations in 
the form of loops of individual fractional-flux vortices are 
responsible for the critical properties of the system in the 
absence of an external field. 

Note that since each individual amplitude is frozen, this 
model will be different from the case where only the sums 
of the squares of the amplitudes are frozen?^ . The latter 
is usually referred to as the A^-component scalar QED 
{NSQED)^^, or the CP^-^ model^S. (As far as criti- 
cal properties are concerned, the NSQED model and the 
(^pAf-i j-QQ(jgi iiave been shown to belong to the same 
universality class^S). We strongly emphasize that we 
must distinguish our model from NSQED and CP(^-i\ 
and will consequently be referring to it as the A^-flavor 
London superconductor (NLS) model. The NLS is in fact 
the natural model to consider for the physical systems 
mentioned in the introduction, in particular pertaining 
to the superconducting mixtures of metallic phases of 
light atoms. As we shall see, the NLS model has physics 
which sets it distinctly apart from the NSQED and the 
QpJV-i jnodels, and it does not have critical properties 
in the same universality class as they do. This becomes 
particularly apparent in the large- A'^ limit, as we shall see 
in section III Dl 



Before we proceed further, it is useful to give another 
form of the action. For brevity we introduce the bare 
phase stiffness of the matter field with flavor index a 
defined as l^/'^")^ = |*[,")|2/M("). Then Eq. ((TJ may 
be rewritten in terms of one charged and A — 1 neutral 
modes as follows (details of this are found in Appendix 
EI). We have 5 = / d^rC, with 



1 ^ 2 

+ 4^ E IV'("fl^^^^P(v(0(")-^^(^))) , (4) 



a, 13=1 



where 



N 



)|2 



(5) 



a=l 



The first term in Eq. Q represents the charged mode 
coupling to the gauge field A, and the remaining terms 
are the A — 1 neutral modes which do not couple to 
A. This means that they have gauge charge equal to 
zero. We will come back to this in Section Hill This form 
Eq. Q will be useful later when we discuss finite field 
effects in section lVIIII We also stress that in the above 
expression should not be confused with defined in 
Eq. (HJ. 

Counting degrees of freedom in Eq. |0J requires care. 
The case A = 1 yields the well known answer that a 
phase variable (which is not a gauge invariant quantity) 
is higgsed into a massive vector field by coupling to the 
vector potential. In the case N = 2, the situation is 
different in the sense that one can form a gauge invari- 
ant quantity by subtracting phase gradients. Thus the 
U{1) X U{1) system may be viewed as possessing i) a lo- 
cal C/(l) gauge symmetry associated with the phase sum 
which is coupled to the vector potential and thus yields 
a massive vector field, and ii) a global C/(l) symmetry 
which is associated with a phase difference where there 
is no coupling to the vector potential. These charged 
and neutral modes are naturally described by the first 
and third terms in Eq. Q), respectively. For N — 3, 
the situation is principally different from both the A = 1 
and A = 2 cases. That is, in Eq. Q for A = 3, we find 
one term describing the charged mode (the first term) 
and three terms describing gauge-invariant neutral phase 
combinations. 

The two neutral modes in Eq. in the A = 3 case, 
cannot be properly described by only two terms, for topo- 
logical reasons. A vortex excitation produces a zero in 
the order parameter space, thus making the superconduc- 
tor multiply connected. A vortex with a non-trivial phase 
winding in any of the three components would result in 
non-trivial contributions to two of three phase-difference 
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terms in Eq. (0J . Hence, for N = 3 an elementary vortex 
i.e. with nontrivial winding only in one of the phases ex- 
cites two neutral modes. In general, when all I-!/)'"-*! differ, 
the bare phase stiffnesses of two neutral modes excited by 
each of the three possible elementary vortices, are differ- 
ent. Thus, the neutral modes in the system are described 
by three phase-difference terms in Eq. Q). These three 
terms are not independent when the condition of single- 
valuedness of each of the N order parameter components 
is enforced, namely that individual phases may change 
only by integer multiples of 2tt around zeroes of the or- 
der parameters. 

Using Eq. (0J as opposed to Eq. I^^l, has advan- 
tages, because the neutral and charged modes are explic- 
itly identified. This facilitates a discussion of the critical 
properties of the iV-flavor system. Moreover, Eq. will 
allow us to identify various states of partially broken sym- 
metry which emerge if an iV-flavor system is subjected 
to external magnetic field2fl,. We will come back to these 
points in detail in Sections IVIII and IVIIII 



D. Vortex representation 

In the following, we transform the model Eq. lO into 
a theory of interacting vortex loops of different flavors. 
The procedure is described in detail in Appendix 
The kinetic energy terms are linearized by introducing 
N auxiliary fields v(")(r). Applying the Poisson sum- 
mation formula and integrating over n^"^ (r) constrains 
the fields v("-'(r) to take only integer values v(")(r). In- 
tegration over all ^("^(r) produces the local constraints 
A • v^") (r) = 0, which are fulfilled by replacing v^"' (r) 
with A X h(")(r) where W-°'\r) are integer- valued fields. 
By applying the Poisson summation once more and sum- 
ming over all h("^(r), the fields h(")(r) take continu- 
ous values h'^"^(r) and the integer- valued vortex fields 
m*^")(r) are introduced. We recognize h(")(r) as the dual 
gauge fields of the theory. To preserve the gauge sym- 
metry of h'^")(r) each vortex field of flavor index a is 
constrained by the condition 



A • m(")(r) = 0. 



(8) 



C. The Villain approximation 

The theory Eq. is discretized on a d = 3 dimen- 
sional cubic lattice with spacing a ~ \ and size , and 
in the phase only approximation the action reads 



N 3 

-13 I V'^"^ P cos(A^6'(") (r) - e.4^(r)) 



a=l 



+ |(V X A(r))2 



(6) 



Here, we have included the inverse temperature coupling 
(3 = 1/T. The symbol A^ denotes the lattice difference 
operator in direction /i in Euclidean space and the po- 
sition vector r runs over all points on the lattice. The 
partition function in the Villain approximation is 



poo N .jr N 

= r,= l nil) 



N 



Y ^^^(A0(") -eA + 27rn("))2 



(7) 



-^(AxAf 



where n'"^ (r) are integer vector fields ensuring 27r period- 
icity, and the lattice position index vector r is suppressed. 
Here, we stress the importance of keeping track of the 27r 
periodicity of the individual phases. For = 1 it has 
been shown that thermal fluctuations in this model ex- 
cite topological defects in form of closed vortex loops. At 
the critical temperature the system undergoes a vortex 
loop proliferation phase transition24i2Si2&. 



Hence, the vortex fields form closed loops. At this stage, 
the action reads 



'y. (A X h("))^ 

A. 2/3|V'(")|2 

TV 



ieA • ^E^ ^ ^^"'^ 



27ri^m(")-h(") + ^(Ax A)2 



(9) 



where the vortex fields m(")(r) are constrained by Eq. 
(HI). We integrate out the gauge field A(r) and get a 
theory in the dual gauge fields h("^(r) and the vortex 
fields m(")(r)ia 



^ = E 



27ri y m(") . h(") -I- y 




(10) 



This generalizes to arbitrary N the results of PeskinS, 
and Thomas and Stoned. In Appendix IeI we generalize 
this result even further by including inter-flavor Joseph- 
son coupling. 

When N > 2 there is an important difference from 
the N — 1 case, which gives rise to entirely novel 
physics. Note how it is the algebraic sum of the dual 
photon fields in Eq. (|10|l that is massive. This differs 
from the case A'^ = 1, where e produces one massive 
dual photon with bare mass e^/2, and the model de- 
scribes a vortex field m(r) interacting through a mas- 
sive dual vector field h(r). However, when N > 2, 
since A • m("^(r) = 0, a gauge transformation h'^")(r) — > 
h(")(r) = h(")(r)-|-A5("'(r) for a G [1,...,N] leaves the 
action in Eq. H1U|) invariant if one of the gauge fields. 
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say h^''^ (r) compensates the sum in the last term in the 
action with Ag('')(r) = - ^7^^ ^S^^Hi")- Thus, even m 
the presence of a gauge charge e, such that the direct 
model is a gauge theory, the dual description is such that 
the individual dual photon fields are also gauge fields. 

Integrating out the dual gauge fields we get a gen- 
eralized theory of vortex fields of N flavors interacting 
through the potential D^°'''^\r) 

N 

Sy ^ TT^ m(")(r)L>("''')(r - r')m('')(r'), 

r,r' aji 

where Sx,y is the Kronecker-delta, and the discrete 
Fourier transform of the vortex interaction potential is 
D("''')(q), given hyi^ 

where A^") = and ^''^ is given by Eq. 

Here, the bare mass mo is the inverse bare 
screening length given by ttiq = e^^*^, and |QqP = 
^j^^-^(2sin((7''/2))^ is the Fourier representation of the 
lattice Laplace operator, where = 2Trn^/L with e 
[1,...,L]. Note that J2a -^^"^ = 1- Note also that when 
= 0, the interaction matrix reduces to 

5("'")(q) = 2/3|V(")py^. (13) 

IQql 

This means that when there is no charge coupling the 
matter fields to a fluctuating gauge field, there is no in- 
teraction between vortices of different flavors. This sim- 
ple case corresponds to Eq. ^ representing a system 
of N decoupled 3Dxy models. Also note that for vor- 
tices of different flavors, t] ^ a, when e 7^ 0, the in- 
teraction matrix tends to vanish when the inter-vortex 
distance is much smaller then the effective penetration 
length A = l/niQ. It follows from the fact that when the 
inter-vortex distance is much smaller than A, the vortices 
interact as if A does not screen, i.e. as if A does not fluc- 
tuate. In this case, it is clear that the action we describe 
is simply that of N decoupled 3Dxy models, i.e. inter- 
fiavor interactions vanish, cf. Eq. (|13|) . For instance, for 
the case N = 2, there will be no interactions between 
vortices of condensate ^q^-* and vortices of the conden- 

('2) 

sate ^'o unless we allow the gauge field to fluctuate. In 
the extreme type-II limit where A — s- cx) only intra-flavor 
interactions between vortices will exist (see also Ref. Isifl . 

The first term of the vortex interaction potential Eq. 
(|12|l is a Yukawa screened potential, while the second 
term mediates long range Coulomb interactions between 
vortex fields. If = 1 the latter cancels out exactly and 
we are left with the well studied vortex theory of the GL 
model which has a charged fixed point for e 7^ 0,2M£. For 



N > 2 we find a theory of vortex loops of A'' flavors in- 
teracting through long range Coulomb with an additive 
screened part. If the number of species A^ grows to in- 
finity and ^'^ — !■ 00, the vortex interaction receives the 
dominant contribution from a diagonal unscreened N xN 
Coulomb matrix. But there are physical situations where 
off-diagonal interactions play an important role even in 
the large- A limit (to be discussed below). One can also 
observe from Eq. (PJ that in the A^ — > 00 limit when all 
components have similar stiffness the magnetic flux en- 
closed by elementary vortices also tends to zero. Thus, 
for N ^ 00 the physics of the model is governed by neu- 
tral modes only. 

The energy density of one straight vortex line of flavor 
a in a distance r larger than the effective penetration 
depth A is found by integrating along the line using the 
last term in the potential Eq. (|12|l only^^. This produces 
an energy term of the form D{r) ~ ln(|r|), and shows 
that such a vortex has logarithmically divergent energy. 

The large- A'^ limit of the NLS serves to illustrate how 
different the physics is from the large-A^ limit of the 
NSQED model and the CP^-^ modeiSLS^. In the large- 
A^ expansion of the NSQED model, only one charged 
fixed point is found (which is infrared stable provided 
2 A > 365), with critical exponent = 1 + 48/A + . . . 
in D = 32i. This is consistent with the results found in 
the large- A^ limit of the CP'^~^) model^^. The origin of 
the difference between these results and the results we 
find for the NLS model is easily traced to the following 
fact. The treatment of the NSQED model in Ref. 27 is 
strictly speaking correct only in the case of type-I super- 
conductivity, since they find that for physical values of 
A^, only a first order phase transition from a supercon- 
ductor to a normal metal takes place (no infrared stable 
fixed point is found for physical values of A^). This is 
correct only for values of the Ginzburg-Landau parame- 
ter K < 0.8/a/2, as has been shown in recent large-scale 
MC simulationSi22i and in earlier analytical treatments^. 
The transitions discussed below where neutral modes ap- 
pear do not significantly depend on whether the system is 
type-I or type-II. Our results are therefore best thought 
of as generalizations to arbitrary A^ of the problem stud- 
ied many years ago by Dasgupta and Halperin on the 
frozen-amplitude A = 1 lattice superconductor model—. 
It is this fact that in the present model the modulus 
of each component is fixed, along with the precise ab- 
sence of internal Josephson coupling between matter field 
species, that brings out the novel physics we shall de- 
scribe, namely the charge-neutral superfluid modes aris- 
ing out of N charged condensate fields. 



E. Dual field theory 

Starting from Eq. (|10|1 the above vortex system may 
be formulated as a field theory, introducing A^ complex 
matter fields 0^"^ for each vortex species, minimally cou- 
pled to the dual gauge fields h^"^. This generalizes the 
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dual theory for = 1 in Refs. '3C|34 The theory reade- - 
(for a comment on the case of general N, see also bottom 
of page 42, Ref. H) 



r Let — 1 



)|2 



1^ 



(a) 



(14) 



'?)|2 



Here, we have added chemical potential (core-energy) 
terms for the vortices, as well as steric short-range repul- 
sion interactions between vortex elements. In the iV = 1 
case, a RG treatment of the term fgh^ yields 



dlnl 



(15) 



and hence this term scales up, suppressing the dual vector 
field h. The charged theory in d = 2-1-1 therefore dualizes 
into a |(/)|"' theory and vice versai^. Correspondingly, for 
N > 2, Eq. (fT^ suppresses X]q=i but not each 

individual dual gauge field. For the particular case N — 
2, assuming the same to hold, we end up with a gauge 
theory of two complex matter fields coupled minimally 
to one gauge field, which was also precisely the starting 
point. Thus the theory is self-dual for N = 3Si&. 



III. CHARGED AND NEUTRAL VORTEX 
MODES 



In this section, we present a straightforward method 
of identifying charged and neutral vortex modes for the 
model Eq. (Q. Consider first the case N = 2, when the 
action Eq. ifTUIl reads 



27ri 



1 

20 



(Vxh(l))2 (Vxh(2))2 



I — 



|V(2)|2 



2/3 



(16) 



From this we identify the massive linear combination of 
the dual gauge fields h^"), namely H = h^^^ + h^^\ If 
a neutral vortex mode exists in the system, this implies 
the existence also of a gauge field in the problem, which 
we will denote by A. Wc therefore write h'^") as linear 
combinations of H and A as follows 



(17) 



We insert this into Eq. H16|l and demand that cross-terms 
between Ti. and A vanish, thus obtaining the following set 
of equations determining the coefhcients (r'^"\ A^^"^) 

p(l) _^p(2) ^ 

A(i)+A(2) = o! 

r(i)AW/|v^(i)|2-i-r(2)A(2)/|v.(2)|2 ^ 0. (i8) 

Thus, we have F^") = It/i^"' |2/^'2^ .^^here *2 = |^(i)|2 
|t/'(2)|2^ which yields the following expression for the 
gauge field A 



A 



1 1^(2) |2h(l) _ 1^(1) |2h(2) 



A(i) 



\I/2 



(19) 



Since we have three equations and four unknowns, we 
may choose A^^^ freely, and determine it by simplifying 
the prefactor in A to get A^^^ = 17*2 = -A(2), whence 
we have 



^=|^(2)|2i^(l)_|^(l)|2i^(2)_ 



(20) 



Inverting the relations for Ti. and A, we have 

h(i) = (|^/.(l)|2H+y^)/*^ 

h(2) = (|^(2)|2^__^)/^2_ 



(21) 



Inserting this back into Eq. 1)16(1 . collecting terms, and 
redefining the fields H/*^ -> Ti and A/'b^ A, we have 
the action S = Su + Sj, where 



Sh = El27ri7i:-m' 



,(+) 



1 



1 



[(V X nf + mlrL^ 



Sa = E|27rU-m(-) + ^(Vx^)2J., (22) 



where 



1 



m(-) = m(i)-m(2), 

(|^(l)|2 + |^(2)|2) 



,(+) 



2/3w 2/3 

1 ^ (l/|7^(l)|2 + l/|V,(2)|2) 

213a ^ 2/3 



(23) 



and ttiq = 62^1*2. The action in Eq. 1(221) . which is equiva- 
lent to Eq. ((16(1 . therefore describes a vortex mode m*^+^ 
interacting with itself via a screened anti Biot-Savart in- 
teraction mediated by the massive vector field Ti, and 
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the vortex mode m^") interacting with itself via an un- 
screened anti Biot-Savart interaction mediated by the 
gauge field A. Hence, the former vortex mode is charged, 
the latter is neutral. In Appendix^ we present an alter- 
native method of identifying charged and neutral modes 
for general N. 



IV. GAUGE FIELD CORRELATORS 

Gauge field correlation functions are useful objects to 
study when considering the critical properties of gauge 
theories. The main reason is that they provide non-local 
gauge invariant order parameters for the theories, which 
in turn enable reliable determination of critical expo- 
nents, including anomalous scaling dimensions. More- 
over, these correlators explicitly identify the mechanism 
by which the Meissner effect is destroyed in type-II su- 
perconductors: The mass of the gauge field A, and hence 
the Higgs phase ( equivalently the Meissner phase ) is de- 
stroyed by a thermally driven vortex loop proliferation of 
the charged vortex mod^^^^^a^^^^. 

In this section, we study in detail the direct gauge field 
correlation function, as well as various combinations of 
dual gauge field correlation functions, in order to gain 
insights into the nature of the critical points Eq. can 
exhibit. 



A. A-field correlator and Higgs mass 

We first consider the propagator for the gauge field A, 
which provides information about at which of the criti- 
cal points the Higgs phenomenon takes place, and where 
the remaining (neutral) fixed points appear. We present 
compact expressions for the general- case, in later sec- 
tions we present explicit numerical results for the cases 
TV = 2 and iV = 3. 

We compute the correlation function (A(r) • A(0)) in 
terms of vortex correlators in the standard way by start- 
ing from the action Eq. lO, prior to integrating out the 
gauge field A, adding source terms containing currents 
J minimally coupled to A, and performing functional 
derivations with respect to the currents that are subject 
to the constraint V • J = 0, after which the currents are 
set to zero. The details of the computations required to 
compute the A-field correlator are given in Appendix lUl 
The discrete Fourier transform of the gauge field propa- 



gator is 5a (q) 



(Aq • A_q). We find 



2//3 



1 



27r^/3e^ G(+)(q) 
IQqP IQqP + m 



, (24) 



where we have defined the correlation function of the 
charged vortex mode as 



N 

GW(q) = ((5:i^(")rm(«)).(5:iV'('')| 

a=l 



N 

El 

r,=l 



ml 



))(-25) 



Notice in Eq. if^ . that the A-field correlator 
is only affected by the gauge-charged vortex mode 



IV' via the couphng constant cx e^. 

Eq. (|24|l is useful in MC simulations, in conjunction 
with scaling forms to be presented below, for extracting 
the gauge field mass and the anomalous scaling dimension 
of the gauge field. The correlation length f a that appears 
in a scaling Ansatz for the A-field correlator 



eA(x) 



1 



Q± 



(26) 



is related to the mass of the gauge field via toa = Ca^ • 
Here, ryA is the anomalous scaling dimension of the gauge 
field A. Consequently, the gauge field propagator Eq. 
(Pl|l has the general structure^ 



eA(q) 



(27) 



|Qq|2 + SA(q)' 
where, close to the critical point 

EA(q)-mi + qqr''A+o(|q|5), (28) 

C is a constant and 5 > 2 ~ tja- By taking the q — > 
limit of the Eqs. H27|l and H28|l we may extract the gauge 
mass from MC simulations. From the relation B = A x 
A the gauge mass is identified as the inverse magnetic 
penetration depth A. The masses of dual gauge fields are 
defined in a similar fashion. 

Let us make a remark concerning how a charged fixed 
point (77A — 1) could be distinguished from a neutral 
fixed point (t^a — 0) by gauge mass measurements. The 
magnetic penetration length is related to the supercon- 
ducting coherence length ^ viaS^t^ 



X- 



(29) 



where v is the critical exponent of the coherence length 
in the superconductor, i.e. v — 0.67155(3)'*^, and d is 
dimensionality. Therefore, we see that when 77A ~ 0, we 
hav e^^i?^ 



while when r/A = 1, we have 



(30) 



(31) 



Hence, the gauge mass ttia = A^^ plotted as a function 
of temperature in the critical regime should for t^a = 1 
give a curve with positive curvature, while for 7/a = it 
should give a curve with negative curvature. 

The compact expression Eq. H24|l is valid for arbitrary 
number of matter field flavors N, and generalizes the 
expression obtained in22i. Note that if ~ 0, we have 
trivially that Eq. (|24|l reduces to 



2/l3_ 

IQqP 



(32) 
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which is always massless. In Sections IVII and IVTl we will 
use large-scale MC simulations to study in detail the case 
N = 2 and = 3, respectively. The main feature of Eq. 
(|24|l is that at low temperatures, we may in the very 
simplest approximation entirely ignore the vortex cor- 
relation function G'^+-'(q) such that t/A(q) is obviously 
massive with photon mass given by the bare mass toq of 
the problem. Actually, in the low-temperature regime, 
we have G^^' (q) ~ q^ which in the long-wavelength limit 
exactly cancels the factor l/|Qqp, rendering the propa- 
gator massive. 

However, at the superconducting critical temperature, 
vortex loops proliferate^SiSSi^SiSLi^ resulting in vortex 
condensation and hence limq^o G^"*"-* (q) ~ const. Now, 
the term inside the brackets in Eq. (|24(l will diverge, 
dominating the behavior of the A-field correlator, such 
that t/A(q) ^ l/l'^- Thus, the Higgs mass is destroyed. 
Note that the amplitudes of the matter fields play no role 



in this, since they are entirely frozen in the present Lon- 
don approximation. It is the condensation of topological 
defects of the matter fields, i.e. vortex loops, that are 
responsible for bringing the Higgs mass to zero, not the 
vanishing of the amplitudes^'' . Therefore, we may view 
the divergence of the penetration length (the correlation 
length in the A-field propagator), as a manifestation of 
the vortex loop blowout in the system. Vortex loops have 
dual counterparts in the current loops of the matter fields 
^'[("■'(r) in Eq. Conversely therefore, we may also 

view the Higgs-mass, i.e. the Meissner effect in the super- 
conductor, as a manifestation of blowout of super-current 
loops upon entering the low-temperature phase. Again, 
the amplitudes of the matter fields ^'g"' (r) play no spe- 
cial role here, other than that they have to be non-zero 
across the Higgs transitio n^^i?^i??i?7i"^9 . 



B. Dual gauge field correlators 



The details of the computations required for finding the dual gauge field correlation functions in terms of vortex 
fields are found in Appendix^ We find the following "Dysons's equation" for the gauge field correlator 

(h(j") • hL^^) = 5("''3)(q) - 7r25("''')(q)5('3''')(q)(m(j'') • mL",^), (33) 



where we have used the fact that the trace of the trans- 
verse projection operator is given by Tr [P!^''] = 2, the 
matrix elements D^"'^\q) are defined in Eq. H12|l . and a 
summation over the indices (j], k) G [1, . . . , iV] is under- 
stood. These results are valid for all N. 

To obtain more explicit expressions, we will work out 
in detail what we obtain for iV = 2. As we have seen 
above, in this case it is natural to use Eq. I|33|) to form 
correlation functions of the combination h^^^ -|~ h'^^^ We 
will, for completeness also consider the combination and 
h(i) - h(2) and IV-^^^ph^i) - IV-^i^f h(2). We also use the 
fact that the interaction matrix D'^°'-'^\q) is symmetric, 
and introduce the definitions 



hW = h(i)±h(2) 

q — q q 



= 5(i^i)(q)±5(i'2)(q) 
= 5(2^2)(q)±5(l'2)(q). 



(34) 



It is enlightening at this stage to introduce the expres- 
sions for (q), as follows 



J(l-l)(q)^^ 
2/3|V'(i)p|i/'(2)|2 

i3(2-2)(q)^2 

i3(l'2)(q)^2 1 ^ 1 



IQqP |^(2)|2|Q^|2 
1 1^(2) |2 1 



IQql^ 



IQc 



(35) 



where *2 = |^(i)|2 _|_ |^(2)|2^ Using Eqs. (ESJ in Eqs. 
(PH), we find 



,(+) 



2/?|V'W|2 



IQqP+m2' 

2/3|V;(2)|2 

|QqP+m2' 



(36) 



and a*- •* and b^ •* given by 



l(-)^2 



2/3|V'(l)|2|?^(2)|2 

b(-)^2 
2/3|?/'(l)|2|V(2)|2 



IQqP 

2 

10^ 



|^(1)|2/|^(2)|2 

IQqP+"l0 

|^(2)|2/|^,(1)|2 



IQc 



-,(37) 



where 



= e2*2 



Notice how the unscreened part 



of the interactions cancel out in (aW,6(+)) but not m 
(a(-),5(-)). This is the origin of the qualitatively dif- 



ferent behavior we will find for the hq^-* and hq ' cor- 
relators. Notice also how the expressions simplify when 
[■^(i)^ = |'(/;(2)|2^ when the screened part of the inter- 
actions appearing in a^^\b^^'> vanish, such that a*-^' = 

We may now write the correlation functions of the two 
relevant linear combinations of dual gauge fields as fol- 
lows 
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I 



(38) 



Using Eqs. (|^ . and we find the surprisingly com- 
pact expression, valid for all N 



2/3 V'^ 



IQq 



(39) 



where we have again introduced G'-"'"-' (q) appearing in Eq. 
(|25|l . In fact, this result could have been written down 
using the known result for the charged case for N — 1^^, 
in combination with Eq. (|22|l . considering the part of 
Eq. H22() only pertaining to the massive vector field H. 
This provides a nice consistency check on the general 
expression for the dual gauge field correlators, as well 
as on the interaction matrix D^°''^^{q). In the low- and 
high-temperature phase, the vortex correlator G*-+-'(q) 
behaves as ^ q^ and ~ c(T), respectively. In either case, 
the dual gauge field correlator Q^"^ (q) is always massive. 

Consider the correlation function of the combination 
of dual gauge fields A = IV'^^^ph^i) - IV'^^^ph^^) which 
couples to the gauge-neutral vortex mode in Eq. (|22|l . 
In principle we may follow the routes used in the above 
calculations, but by now we realize that a quick way of 
obtaining the results is to use Eq. (|22|l in combination 
with the known results for the case = 1 in the neutral 
— . We define 



casi 



(40) 



and find immediately, using the results of Ref. |22 along 
with the definitions in Eq. I|23|) 



GAiq) 



where 



G(-)(q) 



213a 

IQqP 



((m« 



1 - 



2^2/3^G(~)(q) 

IQqP 



.(2) 



) ■ (m' 



(1) 
q 



,(2) 



(41) 



)) (42) 



is the correlation function of the gauge-neutral vortex 
mode. 

In the long wave length limit the behavior of G'-^^(q) 
gives rise to a dual Higgs mechanism. This comes about 
because the G*-~' (q) correlation function is always ~ q^ 
at long wavelengths, but has a non-analytic coefficient in 
front of the q2 term given by the helicity modulus of the 



J 



gauge-neutral mode m'^^^ ~ m'^2)_ This serves to cancel 
the 1/q^ term in the Ga{(1) correlation function exactly. 
This cancellation, originating in the vanishing of the he- 
licity modulus of the gauge-neutral mode, is responsible 
for producing a dual Higgs mass ttia in 5>t(q). Higher 
order terms determine the actual value of the dual Higgs 
mass. Thus, we see that while h'-'^' + h^2) jg always mas- 
sive, iTA^^^ph^-^^ — IV'^^^ph^^^ plays the role of a gauge 
degree of freedom which provides a dual counterpart to 
A in Eq. (Q. This is a manifestation of the self-duality 
of the theory which we have alluded to abov o^i^i^^ . 

Notice that the existence of a dual Meissner effect aris- 
ing out of Eq. (|41() is a substantially more subtle effect 
than the direct Meissner effect coming out of Eq. I|24|l . 
The correlator of the gauge-neutral mode has the prop- 
erty 



GM(q) = Gaq' + G4q^ + 0(q6 



(43) 



for all temperatures, in analogy with the vortex corre- 
lator of the 3Dxy model for the case = 1. It is the 
non-analytic behavior of the coefficient C2 , involving the 
helicity modulus of the gauge-neutral mode, which is re- 
sponsible for producing a dual Higgs mass as the gauge- 
neutral mode proliferates. To obtain a dual Meissner 
effect, a subtle cancellation is required, namely that at 
some critical temperature Td, we must have 



27r^/3|^Wn^(^)pG2(r.i) _ 

|l/.(l)|2 + |V>(2)|2 



(44) 



where we have used the expression for from Eq. H23|l . 
It is important to note that while the actual value of the 
dual Higgs mass is influenced by the higher order terms 
in Eq. (|43|l . the criterion for obtaining a dual Higgs phe- 
nomenon is only determined by the cancellation among 
the terms of order 1/q^ terms in Eq. (|41|l . This differs 
from the mechanism that destroys the Higgs mass in the 
A correlator, since there no such subtle cancellations are 
required, it suffices that the correlator G'-'*'-'(q) changes 
behavior from a constant to ^ q2 in the long-wavelength 
limit. 

We finally consider the correlation function of h'^'. 
Applying the results from Eq. H38() . wc findi^ 



)A(2)*2g(-)(q) 2^2^(;^(l) „ A(2))G(")(q) 



IQqP 



|Qq|2-|-m2 



(AW-A(2))2t;W(q), (45) 



where A'-"^ = |^/;(") |2/xl'2^ v]/2 = |i/)(^)|2 + |^/)(2)|2^ and the mixed gauge-neutral and gauge-charged vortex field 
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correlator is given by 



UK 



(1) 



.(2) 



2 



).(}^|V'(")|W"<^)). (46) 



2 (a)\ 



Note that for the case iV = 1, such that either A'-^-' or 
A*^2) vanishes, then the remaining A('') = 1, only the last 
term in Eq. H45|l survives, and Gi^ "* (q) correctly reduces 
to G^\<i) in Eq. H39|l . In the long wave length limit, 
it is the second term in the curly brackets in Eq. 145|) 
that dominates, giving rise to a dual Higgs mechanism. 
Notice again how it is the vortex correlator G^~^ (q) which 
determines the fate of the massless dual gauge field h^^^ — 
h(2), just like in Eq. 141|l . This is particularly evident 
for the case IV'*'"'^''! — I'tp^^M, when Eq. H45I) reduces to 



'(q) 



4/3|V>(i) 



1- ^2^1^(1) P 



IQqP 



(47) 



This correlator for N ^ 2, e ^ has precisely the same 
form as the dual gauge field correlator for the case N = 

1, e = 0, which exhibits a dual Higgs phenomenonS^. 
Substituting A*") = jV'^^^lV*^ in Eq. (gSJ, we see 

that the criterion for destroying the dual Higgs mass is 
precisely the same as the criterion we arrived at in Eq. 
(I44|l . Thus, whether we compute the correlator in Eq. 
(I45|l or that in Eq. H41|l to establish the existence of a 
dual Higgs phase does not matter. Furthermore, for N — 

2, e 7^ 0, m'^^^ — m(2) behaves as vortices for = 1, e = 0, 
i.e. it is a superfluid mode arising out of superconducting 
condensates. A nonzero mji for the dual gauge field A is 
produced by disordering 0*^^^ at a critical temperature Td 
while a nonzero ttia for the gauge field A is destroyed by 
disordering 0'^2) ^ critical temperature Tc2. 



V. MONTE CARLO SIMULATIONS, iV = 2 

Since the bare interaction between vortices is domi- 
nated at long distances by an unscreened part, it is of 
interest to study the character of the phase transition 
associated with the generation of a Higgs mass for the 
gauge field A. For the = 1 case, it is known that the 
vortex tangle of the 3Dxy model is incompressible and 
the dual theory is a gauge theory such that (</>) ^ is 
prohibited. For the charged case, the vortex tangle is 
compressible, the dual theory only has global symmetry, 
and hence vortex condensation and {(f)) 7^ is possible. 
The introduction of charge destabilizes the 3Dxy fixed 
point. 

To investigate what happens for the case N — 2, MC 
simulations have been carried out for the action Eq. I|ll|) 
on a three dimensional lattice of size L x L x L ioi two 
different cases. In the first case we simulate with unequal 
bare stiffnesses IV-'^^P = 1/2 and IV'^^^P = 1, = 1/4 
and ttIq = 3/8. The bare stiffnesses have been chosen 
to have well-separated bare energy scales associated with 
the twist of the two types of phases. In the second case we 



use equal phase stiffnesses {i/j 



(1)|2 



2)12 



1/4 



and rriQ — 1/2. The values for mo have been chosen such 
that they are of order the lattice spacing in the problem 
to avoid difficult finite-size effects. One MC update con- 
sists of inserting a unitary vortex loop of random direc- 
tion and species according to the Metropolis algorithm. 

To calculate the critical exponents a and v we per- 
formed finite size scaling (FSS) analysis with bootstrap 
error estimates of the third moment of the actionii 
M3 = {{Sv - {Sv)f)/L^ where Sy is given in Eq. 
The peak to peak value of this quantity scales with sys- 
tem size L as whereas the width between the 
peaks scales as L^^^'^. The advantage of this is that 
asymptotically correct behavior is reached for practical 
system sizes. 

To characterize the phase transitions further, we con- 
sider the correlation functions given in Eqs. H24() . H25|l . 
and l|39|) . In the Higgs phase the gauge field mass 
scales according to the Ansata^ given by Eqs. (|77|l and 



f3 



+ C\q\'-^-+0{\qn (48) 



with a corresponding Ansatz for ^^^■'(q)- The masses 

of A and X]^=i therefore defined through the 

q ^ limit of the respective Ansatze 



lim — -T — 

1™ —7T\ ■ 



(49) 



The gauge field masses arc found by measuring vortex 
correlators followed by a fit for small q to their respective 
Ansatze. 

We briefly review the = 1 GL-model. The dual field 
theory of the neutral fixed point is a charged theory de- 
scribing an incompressible vortex tangle^S. The leading 
behavior of the vortex correlator G''^+^(q) ~ (mq-m_q) 

:-22 



lim G(+)(q) 



[l-G2(T)]q2 
q2-C3(r)|q|2+''i. 
l2+C4(r)q4 



T<T, 

T = T, (50) 

T>Te, 



where r]h is the anomalous scaling dimension of the dual 
gauge field h. For T < Tc the mass of the dual gauge field 
given by Eqs. and (gHJ (with = 1 and e = 0) is 
zero, however for T > the q^^ terms in Eq. (|39|l can- 
cel out exactly and the mass rrih attains an expectation 
value. At the charged fixed point of the GL model, the 
effective field theory of the vortices is a neutral theory. 
The vortex tangle is compressible with a scaling Ansatz 
for the vortex correlator 



lim G^^ 



'(q) 



q' ■,T<T, 

|2-,,A ;T = rc 



|qr 

c(T) 



(51) 
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where c(T) is a nonzero constant. Consequently, from 
Eqs. (Egi, and 63) (with = 1 and e ^ 0), the 

mass mA drops to zero at Tc, and the mass of the dual 
vector field nih is finite for all temperatures and has a 
kink at Tc^^. Renormalization group arguments yield 
T]A = 4 — d where d is the dimensionalitjiSL^Siii^ which 
has recently been verified numericalljiSSiiS. 

A. Critical exponents a and u, IV"'^'! < \ip^'^^\ 

We observe two anomalies in the specific heat at Tci 
and Tc2 where Td < Tc2- We find Td and Tc2 from scal- 
ing of the second moment of the action {{Sy — (Sy))'^) / 
to be Tci = 1.4(6) and = 2.7(8). The Mg FSS 
plots for system sizes L = 4, 6, 8, 10, 12, 14, 16, 20, 24 are 
shown in Fig. ^ From the scaling we conclude that 
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FIG. 2: G<+^(q) for iV = 2 and L = 32, plotted for temper- 
atures T = 2.86 > Tc2, T = 2.76 ~ rc2, and T = 2.63 < rc2, 
liniq^o G'"'"' (q) ~ c(r), ~ jqj, and ~ q^, respectively. The 
q —> behavior of the correlator matches precisely the signa- 
ture of a changed fixed point given in Eq. I|51|l . 
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FIG. 1: The FSS of the peak to peak value of the third mo- 
ment AMs labeled (□) and (+) for Td and T^2, respectively. 
The scaling of the width between the peaks A/3 is labeled (A) 
and ( X ) for Td and Tc2 respectively. The lines are power law 
fits to the data for L > 6 used to extract a and v. 

both anomalies are in fact critical points, and we ob- 
tain a = -0.02 ± 0.02 and = 0.67 ± 0.01 for Td and 
a = -0.03 ± 0.02 and ly = 0.67 ± 0.01 for Tc2. These 
values are consistent with those of the 3Dxy and the in- 
verted SDxy universality classes found with high precision 
in Refs. .44.45»4& 



B. Vortex correlator, Higgs mass, and anomalous 
scaling dimension, < l'/''^'! 

The vortex correlators for the N — 2 case are sampled 
in real space and G'"'"-' (q) given in Eq. (|25|l is found by 
a discrete Fourier transformation. At the lower transi- 
tion Td the leading behavior is G''+-'(q) ~ on both 
sides of the transition. Consequently, due to Eqs. (|24|l . 
(|39|l . and (|49|) . to a and to^h are finite in this regime. 
This shows that the vortex tangle is incompressible and 
that the anomalous scaling dimension t^a = 0, which 
corresponds to a neutral fixed point. Fig. [21 shows the 
correlator G''^+-'(q) around Tc2- Below Tc2 the dominant 
behavior is G^~^''{q) ~ q^ whereas G'--^^{q) ~ c(T) above 



the transition. At the critical point G*-+-'(q) ^ |q|, in- 
dicating 77a = 1- Accordingly to a is finite below the 
transition and zero for T > Tc2- 

For each couphng we fit GA{<i)^^ for |Qq| < 0.9 using 
system sizes L = 8, 12, 20, 32 to Eq. (gHJ- The results 
for TOa, and msh which is found in a similar fashion, are 
given in Fig. |3 The system exhibits Higgs mechanism 
when TOa drops to zero at Tc2 with an anomaly in TO^h 
due to vortex condensation. Furthermore toa has a kink 
at Td due to ordering of the phase difference O^^'' — 0^^-' 
with the phase stiffness p ji/'^^) |V(2|V'^^^ P+2|V'^^^ 
confirm Eq. 6J^^. The anomalies in toa and TOsh co- 




FIG. 3: The mass ttia (♦) and 1 — mo/msh (+) found from 
Eqs. ^ and Two non-analyticities can be seen in 

niA at Tci and rc2, corresponding a neutral fixed point and 
a charged Higgs fixed point, respectively. An abrupt increase 
in msh due to vortex condensation is located at rc2. 

incide precisely with Tc2 and Td. Note also how toeh 
changes abruptly at Tc2- This is due to a sudden change 
in screening by h^"-*, giving an abrupt increase in 

TOsh. This is consistent with the flow equation Eq. I|15|l . 
Note that the mass of the algebraic sum of the dual fields 
appears in Eq. H10|) after integrating out the gauge field 
A. 



We may understand the transitions as follows. Above 
Tc2, A is massless, giving a compressible vortex tangle 
which accesses configurational entropy better than an in- 
compressible one. Below Tc2, A is massive and merely 
renormalizes |^'|'' terms in Eq. The theory is effec- 
tively a I^E'I'* theory in this regime. Thus, the remaining 
proliferated vortex species originating in the matter fields 
with lower bare stiffnesses form vortex tangles as if they 
originated in a neutral superfluid. For the general N case, 
a Higgs mass is generated at the highest critical tempera- 
ture, after which A renormalizes the term, such that 
the Higgs fixed point is followed by — 1 neutral fixed 
points as the temperature is lowered. 

The picture that emerges from the above discussion 
of the gauge field and the dual gauge field correlators 
is the following. Below Td there is one massless "pho- 
ton" , namely j^^^^^ph^^^ — l'!/;^-^^ ph(^\ while A is massive. 
Above Tci and below rc2, both jV'^^^ph^i) - ji/'^i) ph^^) 
and A are massive, while above Tc2, ph*^^^ — 

|^(i)|2]^(2) jg massive and A is massless. 

C. Critical exponents a and u, \tl>^^'^ \ = |^''^'| 

A special case is obviously presented by the case 
1^(1) I ^ 1^(2) I gjj^(,g ^j^gj^ j,^^ ^ j.^^ ^ g^j^j Yva-ve 

a transition directly from a low-temperature phase with 
one massless dual gauge field ji/j^^^ph^^^ — ji/j^^^ph^^^ = 
1^(1) |2(']^(i) _ jjCs)-) high-temperature phase with one 
massless direct gauge field A. Thi s is the remarkable 
self-duality observed in Refs. I5l l6ll9i 

The second moment of the action with — 
|'0(2)|2 = 1^ = 1/4 and TOq = 1/2 exhibits one anomaly 
at Tc — 2.7(8). Scaling plots of the third moment of the 
action are shown in Fig. ^ FSS yields a = 0.03±0.04 and 
V = 0.60 ± 0.02. The numerical value for v is in agree- 
ment with the value found in Ref. HI, = 0.60 ± 0.05. 
Note that our result for a and v is not in agreement with 
hyper scaling. 
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FIG. 4: The FSS of tlie peak to pealc value of the third 
moment AM3 labeled for Tc and for = The 

scaling of the width between the peaks A/3 is labeled (o). The 
lines are power law fits to the data for L > 8 used to extract 
a and i/. 
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D. Vortex correlator and Higgs mass, \ip^^''\ = IV'^'^^I 

The mass of the gauge field ttja was found by fitting 
^A(q)^^ data from system sizes L — 8, 12, 20, 32 to Eq. 
(j48|l . The mass m^h was found similarly. The results are 
presented in Fig. El 
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FIG. 5: The mass jtia (♦) and 1 — mo/mEh (+) found from 
Eqs. ^ and (EHJ, for = \il>^^^\. One non-analyticity 

can be seen in niA at Tc, corresponding to a fixed point which 
is not in the 3Dxy or inverted SDxy universality class. An 
abrupt increase in msh due to vortex condensation is located 
at Tc = 2.7(8). 



E. Discussion 

The result for the exponents a and f at for \^^^^ \ = 
|'(/''^^)| shows that when the SDxy and inverted SDxy crit- 
ical points collapse onto each other, then instead of a 
simple superposition, one gets a new fixed point which 
is in a different universality class. This result is far from 
obvious. Naively one would perhaps have guessed from 
Eq. H22() that for = 2 one has two decoupled vortex 
modes, one neutral mode exhibiting a phase transition 
in the SDxy universality class and one charged mode ex- 
hibiting a phase transition in the inverted SDxy univer- 
sality class. At I = I a naive guess would be that 
one would have two such phase transitions superimposed 
on each other, giving a and v values in the SDxy uni- 
versality class. However, there is a principal distinction 
from the case when IV'^'^-'I ^ IV'*'^''!- In the latter case the 
upper phase transition is always a charged critical point 
because the neutral mode is not developed. Thus at the 
upper transition the interaction of vortices is of short 
range, while at the lower transition there is a prolifera- 
tion of vortices with long range interaction. However, in 
the case l?/'^^-' | = I? then below the single phase tran- 
sition both types of vortices have neutral vorticity along 
with charged vorticity and thus this phase transition can 
not be mapped onto a superposition of a neutral and a 
charged fixed points. 

Also, it is the fact that the system is self-dual at this 
point that invalidates the naive superposition conjecture. 
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since the 3T)xy and inverted 3Dxy phase transitions do 
not describe phase transitions of a self-dual system. Even 
though the value of v appears to be in good agreement 
with the 3D Ising value, we observe that the 3D Ising 
model is not self-dual either, and the new type of criti- 
cal point for = can therefore not be in the 
3D Ising universality class. The origin of the novel ex- 
ponents is therefore essentially topological, showing that 
when the vortex loop blowouts of the neutral and charged 
modes are not well separated, they interact in a non- 
trivial fashion. There will therefore exist a crossover 
regime parametrized by the field It/j^^^P — where 
the exponents a and v change from SDxy values to the 
new values we find here (see Fig. 7 of Ref. 0) • In princi- 
ple, it is possible to compute the relevant crossover expo- 
nents in order to shed further light on this new self-dual 
universality class. 



VI. MONTE CARLO SIMULATIONS, = 3 

In the model Eq. I|ll|l with = 3 vortex flavors we 
expect in general one charged critical point associated 
with the condensation of the charged vortex mode and 
two neutral critical points where neutral vortex modes 
proliferate. To study the phases of this model we have 
performed MC simulations with the action given in Eq. 
(Ill|l with bare phase stiffnesses given in Tab. D We 
have applied the same methods for calculating the critical 
exponents a and v as well as gauge masses as we did for 
the N = 2 case. 

It is useful to give the superfluid modes specifically 
for the = 3 case (see Appendix El for details of the 
derivation for the general- TV case). Using Eq. we 
have for this case 



S = 



1 / IV-^I^ 



*2 I 

|V.(l)p|^(2)|2 
2*2 



2)12 



1^(3) I 



)2 
1 
+ y({V'(")})-h-(Vx A)^ 



(V(^^(l) - ^^(2)))2 + M^W!^(V(0(1) - 0(3)))^ + M!!W!^(V(^?(^) - 0(3)))^ 



2*2 



2*2 



(52) 



Here, we have defined *2 ^ |^(2)|2 _^ |^(3)|2_ ^^le regime of short penetration length, the combination 

of phase gradients which is coupled to the gauge field A can be gauged away at length scales of the order of the 
penetration length A — 1/e*. The remaining gradient terms for the neutral modes are given by 



2*2 



2*2 



2*2 



(53) 



This action could be inferred also directly from Eq. Hll() . For the case iV = 3, we write the action in the vortex 
representation as 



3 3 



q r]—l a—1 

which when written out takes the form 



Sv _^/(EJ^(")pm("').(EjV'('')pm^'^) , |^Wn^(2)p(mW-m(^V(mLV«^-q) 



2^2/3/*2 ^ m^ + ml IQql' 



q 



1^(1) |2 1^(3) |2(^(l) _ ^(3)) . (^W _ ^(3)) |^(2)|2|^(3)|2(^(2) „ ^(3)) . (^(2) _ ^(3)) 
IQqP ^ IQqP 



(55) 



The three last terms in Eq. (|55|l are nothing but the 
vortex representation of Eq. H53() . Notice also how all 
cross-terms between different vortex species cancel out 
for arbitrary bare phase stiffnesses when TOq — 0. 

Thus, for the case TV = 3, we have three phase vari- 
ables yielding three neutral gauge invariant combinations 



of phase differences. This amounts to two true neu- 
tral modes, the remaining degree of freedom is associ- 
ated with the composite charged mode, which absorbs 
A and yields a massive vector field via the Higgs mech- 
anism. If all three bare phase stiffnesses IV'^'^-'Ij 
and |?/;*-3)| are different, this yields one charged inverse 
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SDxy critical point where the Meissner effect sets in, 
and two neutral SDxy critical points at lower tempera- 
tures, all separate. Consider now 1-0'"^'! = < l^/"^^^!- 
Then the charged mode proliferates at the highest criti- 
cal temperature where the Meissner-effect sets in, and the 
two neutral modes proliferate simultaneously at a lower 
temperature. The highest transition is still an inverted 
SDxy transition, the lower one is a neutral SDxy criti- 
cal point. Note how this is dramatically different from 
the case N = 2, when the original neutral SDxy critical 
point was collapsed on top of the inverted SDxy critical 
point, resulting in a new universality class of the phase 
transition, essentially due to the self-duality of the N = 2 
system. It is also evident that collapsing a neutral and a 
charged fixed point is quite different from collapsing two 
neutral fixed points. 

For the case = IV^'^'l < l"'/''^'!, in terms of the 

masses of A and the two dual gauge fields associated 
with the neutral modes, rriA is non-zero below the up- 
per critical temperature, while the two dual gauge fields 
become massive above the lower critical temperature. In 
this case, the degenerate lower critical point is therefore 
a SDxy critical point, while the upper critical point is an 
inverted SDxy critical point. 

A further interesting possibility is to set < 
_ 1^(3) 1^ Consider the masses of A and the two 
dual gauge fields associated with the neutral mode in 
this case. At the lower critical temperature, one neu- 
tral vortex mode proliferates in a SDxy transition, gen- 
erating a mass to the dual gauge field (thus breaking 
one dual gauge symmetry). This mode is therefore dual- 
higgsed out of the problem at higher temperatures. The 
gauge field A becomes massive below the upper critical 
temperature, while the dual gauge field associated with 
the remaining neutral mode becomes massive above the 
same upper critical temperature. Hence, the situation 
at the upper critical point corresponds precisely to the 
case N = 2, IV'*-^-'! — IV'^^^'L for which we have already 
seen that a non-SDxy critical point emerges. When all 
bare stiffnesses are equal, IV'^"'"^! = IV'''^''! = IV'^^^li 
three fixed point collapse. We present MC simulations 
for the three cases given in Tab. Q of which the case 
IV'*-^-'! < IV'*'^''! < li'^^M is the most pertinent to mixtures 
of superconducting condensates of for instance hydrogen 
and deuterium, or hydrogen and tritium. 



TABLE I: Phase stifTnesses and bare masses mo for the 

= 3 MC simulations. In all cases the charge e — 1/2. 



Case 


1^(2) |2 




mi) 


1 1/3 


2/3 


4/3 


7/12 


2 1/2 


1/2 


4/3 


7/12 


3 7/9 


7/9 


7/9 


7/12 



A. Critical exponents a and u, < \ip^^^ < \ip^^^\ 

MC simulations are performed for a iV = 3 system 
with bare phase stiffnesses j^'^^^P = 1/3, [V^^^^p = 2/3, 
1^(3) |2 ^ system sizes L = 4,6,8,10,12,14,16. 

We sample the second moment of the action Eq. (|ll|l and 
find three anomalies for temperatures Td, Tc2, and Tcs, 
which from FSS are found to be Td = 0.98, Tc2 = 1.92, 
and Tc3 = 3.63. 

From a FSS analysis of the third moment of the action, 
we have measured the critical exponents a and v. The 
FSS plots are given in Fig. M We find a = -0.03 ± 0.02 
and z/ = 0.65 ± 0.02 for Td, a = -0.02 ± 0.02 and v = 
0.66 ± 0.01 for Tc2, and a = -0.01 ± 0.03 and = 0.69 ± 
0.02 for Tc3. These values are consistent with the values 
for the SDxy and the inverted SDxy universality classes. 




FIG. 6: FSS of the peak to peak value of the third moment 
of action AMs for iV = 3 with {^p'^^^f = 1/3, = 2/3, 

|^(3)|2 ^ labeled (A), (A), and (.), for Td, Tea and T.s, 
respectively. The scaling of the width between the peaks A/3 
(o), (■), and labeled (□), for Td, Tc2, and Tcs respectively. 
The lines are power law fits to the data for L > 6 used to 
extract a and v. 



B. Vortex correlator, Higgs mass, and anomalous 
scaling dimension, Ii/j'^'I < |V''^'| < 1^*^' 

In the Higgs phase, we expect the gauge field correlator 
Ga{(i) in Eq- (|24|) to scale according to the Ansatz Eq. 
(|48ll . For each coupling we fit Ga{(i)~^ from the MC 
simulations for system sizes T = 8, 12, 20 and estimate 
the gauge field mass toa- 

The results for the vortex correlator G^+'(q) in Eq. 
(|25|l and the Higgs mass Eq. H49|) are given in Fig. [7| 
Note how the q-dependence of G^+-'(q) changes when 
the temperature is varied from above to below T^s from 
G'(+)(q) ~ const to G^+^(q) ~ g^, respectively. Note 
also how the q-behavior of the vortex correlator remains 
unchanged when the temperature is varied through Tc2 
and Td, i.e. it remains G(+Hq) q^- This reflects the 
fact that the field A has been higgsed out of the problem 
at Tc3 such that the vortex tangle is incompressible below 
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this temperature. From Eq. (|49() it is therefore clear that 
a Higgs mass is generated at T^s by the estabhshing of a 
charged superconducting mode. Moreover, when the two 
additional neutral superfiuid modes are established at Tc2 
and Tcij this adds to the total superfiuid density and 
hence leads to kinks in the London penetration length 
and thereby ttia- 

Precisely at Tcs, uia. vanishes, and the scaling Ansatz 
given by Eq. H51I) may be used to extract tja.. From Fig. 
□ and G(+)(q) at Tcs, we extract ryA = 1, from which we 
conclude that the critical point at T^s is an inverted SDa;?; 
critical point. Likewise, from the G^"'"^(q) ~ behavior 
at Tci and Tc2 we conclude that these two critical points 
feature 77A = and hence represent 3T)xy critical points. 



C. Critical exponents a and z/, \ip^^^\ = \ip^^^\ < l^'^'l 

MC simulations have been performed for a = 3 
system with bare phase stiffnesses = = 

1/2 and = 4/3 and system sizes L = 

4, 6, 8, 10, 12, 14, 16. By measuring the second moment 
of the action Eq. (|ll|l we find two anomalies for the tem- 
peratures Tci and Tc2, which from ESS are found to be 
Tci = 1.46 and Tc2 = 3.63. From a ESS analysis of the 
third moment of the action we have measured the critical 
exponents a and ly. The ESS plots are given in Fig. |S1 
We find a = -0.03 ± 0.02 and ly = 0.65 ± 0.02 for Td, 
and a = -0.03 ± 0.03 and ly = 0.68 ± 0.02 for Tc2- These 
values are consistent with the values for the SDxy and 
the inverted SDxy universality classes. 
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FIG. 7: Results for the vortex correlator Eq. I25L and the 
Higgs mass Eq. (gHJ for the case = 1/3, = 

2/3, = 4/3. The upper panel shows G^+^(q) as a 

function of |Qq| for seven temperatures starting from above: 
Above and close to Tc3 , above and close to Tc2 , above and close 
to Tci , and below Td . Above Tcs , the vortices are seen to have 
condensed, G^'^^q) ~ const while for all temperatures below 
Tc3 , including above and below Tci and Tc2 , G*-^' (q) ~ for 
small q. The lower panel shows the Higgs mass as a function 
of temperature, showing the onset of Meissner effect at Tcs , 
and the additional anomalies at Td and Tc2 due to the ap- 
pearance of additional neutral modes at these temperatures. 
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FIG. 8: FSS of the peak to peak value of the third moment 
of action AM3 for iV = 3 for = \tp^^^\^ = 1/2 and 

j^^(3)j2 ^ 4^3^ labeled (a), and (•), for Td and Tc2, respec- 
tively. The scaling of the width between the peaks A/3 labeled 
(o), (■), for Td and Tc2 respectively. The lines are power law 
fits to the data for L > 6 used to extract a and v. 



D. Vortex correlator, Higgs mass, and anomalous 



',(3)1 



scaling dimension, \tp^ '\ — \tp^ ' \ < lip 

Like the previous case, we extract the gauge field mass 
by fitting the gauge field correlators for small q to the 
Ansatz Eq. for system sizes L = 8, 12, 20. 

The results for the vortex correlator G'^"'"-'(q) in Eq. 
(|25|1 and the Higgs mass defined in Eq. (|49|l are given in 
Fig. Note how the q-dependencc of G^+-'(q) changes 
when the temperature is varied from above to below 
Tc2 = 3.63 from G(+'(q) - const to G(+)(q) - q^, re- 
spectively. Note also how the q-behavior of the vortex 
correlator remains unchanged when the temperature is 
varied through Td = 1.46, i.e. it remains G*^+^(q) ~ q^. 
This reflects the fact that the field A has been higgsed 
out of the problem at Tc2 = 3.63 such that the vortex tan- 
gle is incompressible below this temperature. From Eq. 
(|49|l it is therefore clear that a Higgs mass is generated 
at Tc2 by the establishing of a charged superconducting 
mode. Moreover, when the two additional neutral su- 
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perfluid modes are established at Tc2 this adds to the 
total superfluid density and hence leads to a kink in the 
London penetration length and thereby toa- 

Precisely at the charged transition Tc2, mA vanishes 
and we find the gauge field correlator has the form 
G(+)(q) - Iqp"''*. From the G(+)(q) data in Fig. Elwe 
extract ?7a — 1 , from which we conclude that the critical 
point at Tc2 is an inverted SDxy critical point. Likewise, 
from the behavior of G^+-'(q) ~ at Td conclude that 
this critical point features t^a = and hence represents 
a 3Dxy critical point. 
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FIG. 9: Results for the vortex correlator Eq. 1125^ . and the 
Higgs mass Eq. ^ for the case It/j'^'P = = 1/2 and 

|.j^(3)|2 _ rpj^g upper panel shows G^'^\q) as a function 
of IQql for five temperatures starting from above: Above and 
close to Tc2, above and close to Td, and below Td. Above 
Tc2, the vortices are seen to have condensed, G^^\q) ~ const 
while close to Tc2, G''^-'(q) ~ |q|. For all temperatures below 
Tc2, including above and below Td, G^'^\q) ~ for small q. 
The lower panel shows the Higgs mass as a function of tem- 
perature, showing the onset of Meissner effect at Tc2, and an 
additional anomaly at Td due to the appearance of additional 
neutral modes at this temperature. 



E. Critical exponents a and 



^{i)| ^ ^ 



(2) 



/,(3) 



MC simulations are performed for a = 3 system with 
equal bare phase stiffnesses = jV'^^-'P = I'i/'^^^p = 

7/9 and system sizes L = 4, 6, 8, 10, 12, 14, 16. From mea- 
surements of the second moment of the action Eq. I|ll|) 



we find one anomaly for temperature the Tc, which from 
FSS is found to be = 2.19. From a FSS analysis of 
the third moment of the action we have measured the 
critical exponents a and v. The FSS plots are given in 
Fig. Uni We find a = 0.02 ± 0.03 and = 0.59 ± 0.02. 
The values appear not to agree with hyper scaling. They 
are not consistent with the 3Dxy universality class. 

The above values for a and j/ are however in agreement 
with those found for the case N — 2, ji/'^^^l = IV'^^^I- We 
observe, based on the numerical results for the two cases 
A^ = 2, IV'^i^l = l^^'^M and A^ = 3, ji/'^^^l = IV'^^^I = IV'^^^I 
compared to the other cases that we have considered, 
that collapsing two neutral critical points in the SDxy 
universality class leads to a single critical point also in 
the SDxy universality class. On the other hand, it ap- 
pears that collapsing A^ — 1 neutral critical points in the 
SDxy universality class and one charged fixed point in the 
inverted SDxy universality class leads to an A^-fold degen- 
erate single critical point in a universality class (which in 
principle depends on A^) which is not that of the SDxy 
or inverted SDxy type. For A^ = 2, we may define the 
universality class as that of a 3D self-dual U{1) x U{1) 
gauge theory, while it is less clear what it is for other 
A^ > 3. 



10 




10' 



FIG. 10: FSS of the peak to peak value of the third moment 
of action AMs for AT = 3 for \tp^'^'>\^ = |i/>'2)p = = 7/9 

labeled (a). The scaling of the width between the peaks A/3 
labeled (o). The lines are power law fits to the data for L > 6 
used to extract a and v. 



F. Vortex correlator, Higgs mass, and anomalous 
scaling dimension, IV'^^'l = i'/'^^'l = l^'^'l 

We extract the gauge field mass toa by fitting the 
gauge field correlators for small q to the Ansatz Eq. 148|) 
for system sizes L = 8, 12, 20, 32. 

The results for the vortex correlator G^+^(q) in Eq. 
(|25|l and the Higgs mass defined in Eq. (|49|l are given in 
Fig. ^2 Note how the q-dependence of G*^+-' (q) changes 
when the temperature is varied from above to below Tc — 
2.20 from G^+Hq) ~ const to G*^+''(q) ^ q^, respectively. 
From Eq. (|49|l it is therefore clear that a Higgs mass is 
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generated at Tc = 2.19 by the establishing of a charged 
superconducting mode. From G*-^^ (q) measurements at 
Tc we find the anomalous scaling dimension to be 77A — 1 • 
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FIG. 11: Results for the vortex correlator Eq. (1251 . and the 
Higgs mass Eq. ^ for the case = IV'^^'P = l^'^'P = 

7/9. The upper panel shows G''''''(q) as a function of |Qq| 
for temperatures above and close to Tc, and below Tc. Above 
Tc, the vortices have condensed, G^^'(q) ~ const. Below 
T, G'^^^(c\) ~ for small q. The lower panel shows the 
Higgs mass as a function of temperature, showing the onset 
of Meissner effect at Tc. 



G. General N 



is a charged fixed point. Proliferation of the vortex 
loops (A6'(i) = 27r, A6'(2) = 0) at Td eliminates the neu- 
tral mode. On the other hand, the composite vortices 
(A^^i) = 27r, Aflf^) = 27r) do not feature neutral vortic- 
ity at any temperature and thus can be mapped onto 
vortices in a A'^ = 1 superconductor with bare phase 
stiffness It/j'-'^'P + j^/j^^^p. A characteristic temperature 
of proliferation of such vortex loops is higher than rc2, 
which excludes the composite vortices from the sector of 
critical fluctuations in the system. The same argument 
applies to the N > 2 case. 
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FIG. 12: (Color online). Phase transitions in the A''-fiavor 
London superconductor with different bare stiffnesses of the 
order parameter components. The green line is the gauge 
field mass ttia- At the highest temperature the system be- 
comes superconducting via a phase transition in the inverted 
3Da;y universality class. At the lower transitions the system 
develops composite neutral superfluid modes in the supercon- 
ducting state via a series of A^ — 1 phase transitions, all in the 
SDxy universality class. 



The critical properties of the A'^-component system are 
governed solely by excitations of vortex loops with frac- 
tional flux. That is, in the N = 2 case, T^i is governed 
by proliferation of the vortex loops with phase wind- 
ings (A6l(i' = 27r,A6'(2) = 0), while Tc2 marks the on- 
set of proliferation of the loops of vortices with wind- 
ings (A0(^) = 0, A(?(^^ = 2n). Remarkably, for general 
N, below the temperature TcAr_i, where Td < ■ ■ ■ < 
TcN-i < TcN, topological excitations with nontrivial 
windings only in one phase has a logarithmically diver- 
gent energjsiSiiS. Moreover, the composite vortex loops 
(A6'(i) = 27r, A6'(2) = 2tt) which in contrast have finite 
energy per unit length, do not play a role as far as critical 
properties are concerned. 

For the case N = 2, the critical point at Tc2 > Td 



Summarizing the previous two sections, the resulting 
schematic phase diagram of the A'^-flavor London super- 
conductor in the absence of external field is presented in 
Fig. El Assuming the bare stiffnesses have been cho- 
sen to have well separated bare energy scales associated 
with the twist of phases of every flavor, we find N dis- 
tinct critical points. At the highest critical temperature, 
the charged vortex mode condenses and the gauge field 
acquires a mass, driving the system into a superconduct- 
ing phase. For lower critical temperatures, neutral vor- 
tex loops condense and the system develops superfluid 
modes. Hence, in zero magnetic field there are A^ — 1 
superfluid modes arising in a superconducting state. 
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VII. N = 2 SYSTEM IN AN EXTERNAL 
MAGNETIC FIELD, LATTICE AND 
SUB-LATTICE MELTING, AND METALLIC 
SUPERFLUIDITY 

We next discuss the situation when the system is sub- 
jected to an external magnetic field. Two important 
aspects of the physics to be described below, are i) 
three dimensionality and ii) a, significant difference in 
the bare stiffnesses of the condensates. As discussed 
recentlji22i2i, when an external magnetic field is applied 
to a three dimensional type-II A^-component supercon- 
ductor, it changes its properties much more dramatically 
than in the ordinary = 1 case. The composite charged 
vortices have finite energy per unit length and couple to 
the magnetic field, and hence are relevant for magnetic 
properties. If the bare stiffnesses of the fields are differ- 
ent, the existence of composite purely charged vortices 
results in a particularly rich phase diagram with several 
novel phases and phase transitions. Note that in the fol- 
lowing two chapters we denote a constituent vortex orig- 
inating in a 27r phase-winding in 9^°'^ a type-a vortex, 
where a G [1, . . . , A^]. 



A. N = 2 system in external field at T = 

In the presence of an external magnetic field, but in 
the absence of thermal fiuctuations, the formation of an 
Abrikosov lattice of non-composite vortices is forbidden 
because these defects have a logarithmically divergent 



energ; 
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cf. discussion following Eq. (|13|l . In a type- 
II A^-component system, the system forms a lattice of 
composite vortices for which A^^"^ = 2tt for every a £ 
[1, . . . , A^]. A schematic picture of the resulting lattice of 
composite vortices in an A'^ = 2 superconductor is shown 
in Fig. 1131 In the discussion below we consider the type- 
II limit, but not extreme type-II since the interaction 
between vortices of different species is depleted at the 
length scales smaller than the penetration length, cf. Eq. 
(|12|l and the discussion following Eq. H13|) . We do not 
discuss effects of this depletion assuming a moderately 
short penetration length scale. 



B. Effects of low-temperature fluctuations on 
field-induced composite vortices 

In this subsection, we will consider the effects of ther- 
mal fiuctuations, and how it affects the Abrikosov vortex 
lattice of composite vortices defined above. 



1. Thermal generation of loop-like splitting of line vortices 

At finite temperature, the A'^-component system sub- 
jected to a magnetic field B will exhibit thermal exci- 
tations in the form of vortex loops with fractional flux 



FIG. 13: (Color online). A type-II, N ^ 2 system at 
zero temperature in external magnetic field forms a lattice 
of composite Abrikosov vortices. A composite vortex may be 
viewed as co-centered type-1 (red) and type-2 (blue) vortices 
(Ae(^) =27r,A6l(2) =27r). 



similar to the B = discussion in the first part of this 
paper. We observe that since the field-induced composite 
vortices are logarithmically boundiSiiS, thermal fiuctua- 
tions will induce a local splitting of composite vortices in 
a configuration of two half-loops connected to a straight 
lin^^tSL as shown in Fig. E| We observe that every 




Two joint 
"half-loops" 




FIG. 14: (Color online). Low-temperature fluctuations in 
the N = 2 system subjected to a magnetic field. Thermal 
fluctuations generate closed loops of composite fractional flux 
vortices and local splitting of field-induced composite vortex 
lines. The type-1 vortices (red) are the vortices of the com- 
ponent with the lowest bare phase stiffness. When these vor- 
tices are viewed as world-lines of bosons, they constitute the 
"lighter" of the vortex species. These "light" vortices fluctu- 
ate more strongly than the "heavier" type-2 vortices (blue). 
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branch of a "split loop" formed on a field-induced vor- 
tex line features neutral as well as charged vorticity. The 
interaction between these two branches is mediated by 
a neutral vortex mode exclusively associated with the 
phase difference 7 — 6^^^ — 9^^"^ . The screened charged 
mode does not contribute to the interaction between the 
two branches. This is implicit in Eq. H12|l as follows. 
The vortex segments of different flavors do not interact 
at short distances much smaller than A = m^^, where 
the charge, or mg appearing in the interaction matrix Eq. 
(|12|l . can be ignored. On such length scales, the screened 
part of the interaction matrix is essentially unscreened, 
and is canceled by the inter-flavor interaction, which is 
unscreened on all length scales. Hence, as also discussed 
in Sectional the intra- vortex interaction is strongly re- 
duced at length scales smaller than A. 

Moreover, in terms of the field 7 = 6*^^^ —9'^^\ two split 
branches of a composite field-induced vortex have oppo- 
site vorticities (A7 = 2tt on one branch and A7 = — 27r 
on another branch). On the other hand, such a loop 
emits two integer flux vortices at its top and bottom, 
which do not feature neutral vorticity. So the process 
of such a thermal local splitting of a field-induced line 
may be mapped onto a thermally generated proliferation 
of closed vortex loops in the artificial phase field 7 as 
those in the neutral N = 1 model in absence of mag- 
netic field'^'^^. Hence, somewhat counterintuitively such 
a splitting transition should be in the 3T)xy universality 
class^liiSi . This transition, being topological in its origin, 
should not be confused with the topological Kosterlitz- 
Thouless transition known to occur in planar systems. 



2. Melting 

Apart from the splitting of composite vortices and gen- 
eration of closed vortex loops, the thermal fluctuations 
will produce one more competing process. That is, the 
lattice of composite vortices can be mapped onto an or- 
dinary vortex lattice in a one-component superconduc- 
tor. Sufficiently strong thermal fluctuations drive a first- 
order melting transition of the field-induced Abrikosov 
latticd^^i^^i^. A counterpart to this effect for the case 
N = 2 when IV''""^'! 7^ IV'^^'I is much more complicated. 
We next consider this process in the regimes of low and 
high magnetic fields, separately. 



C. Sublattice melting in low magnetic fields 

Consider the case of weak magnetic field (much smaller 
than the upper critical magnetic field for which supercon- 
ductivity is essentially destroyed) for the situation where 
\^p^^^ \ IV'^^^I- Introducing a characteristic temperature 
associated with a melting of the type-2 vortex lattice in 

the absence of the condensate ^'q^'', then at sufficiently 
low magnetic field this melting temperature will be much 



higher than the characteristic temperature of thermal de- 
composition of a composite vortex line into two individ- 
ual vortex lines. Thus, the first transition that would 
be encountered upon heating the system, is the thermal 
splitting of field-induced composite vortices into sepa- 
rate type-1 and type-2 vortices. This would be accom- 
panied by a proliferation of closed loops of type-1 vor- 
tices, while the vortices of type-2 will remain arranged in 
a lattice. We will denote this phase transition as sublat- 
tice meltino^Si^. The critical temperature of this phase 
transition is denoted Tslm (see Fig. [TH|) . A schematic 
picture of the sublattice vortex liquid is given in Fig. 
[TKl As discussed above, upon thermal decomposition of 
the composite vortices, the emerging individual vortices 
can be mapped onto positively and negatively electrically 
charged strings which logarithmically interact with each 
other. 

Quite remarkably, the Abrikosov lattice order for the 
component with the highest phase stiffness survives the 
decomposition transition, for the following reason. The 
dominant interaction between individual vortices is the 
long-ranged interaction mediated by neutral vorticity, cf. 
Eq. (|12|l . This permits a mapping of such vortices onto 
positively and negatively charged strings. Upon thermal 
decomposition, the effective long-range Coulomb interac- 
tion mediated by the neutral mode is screened without 
affecting the charged modes. Consider the case when 
\ <C IV'^^''!- Then the stiffness is large enough 

to keep the type-2 vortices arranged in a lattice while the 
stiffness is too weak to constrain type-1 vortices to 

the lattice. Thus, the "light" type-1 vortex lines are in 
their molten phase. This is the physical origin of the 
sublattice melting process. The situation is illustrated in 
Fig. El We emphasize that the existence of the regime 
of sublattice melting follows from the fact that the stiff- 
ness of the neutral mode, which keeps composite vortices 
bound at low temperatures, is always smaller than the 
smallest stiffness of the individual condensates, namely 

- |^(1)|2^.|^(2)|2 < 1^ I • ^^'^) 



D. Composite vortex lattice melting in strong 
magnetic fields 

It is known from the = 1 system that an increase in 
magnetic field suppresses the melting temperature of the 
vortex lattice^. Thus, an important and characteristic 
feature of the phase diagram of the N = 2 system is 
that the composite vortex lattice melting curve should 
at some point cross the decomposition curve. Thus, the 
phase diagram should feature a composite vortex liquid 
phase in the low-temperature, high-magnetic field corner. 
A schematic picture of this phase is given in Fig. E| 
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FIG. 15: (Color online). A vortex liquid of type-1 vortices 
(red) immersed in a background of a type-2 vortex lattice 
(blue) in the N = 2 system in the regime |'(/''^'! <^ i^'^'l- 
This is the type-1 vortex sublattice melting. There is a tem- 
perature region in low magnetic field when "light" vortices 
are decoupled and form a liquid. "Light" vortex loops are 
proliferated, while "heavy" vortices form a lattice immersed 
a hquid of "heavy" vortex loops. Both heavy and light vor- 
tices carry only a fraction of magnetic flux quantum in this 
state. 



(2) 

Ae =2n 



Ae'^=2K 




FIG. 16: (Color online). Liquid of composite vortices in the 
N = 2 model immersed in a liquid of non-proliferated vortex 



loops, 
fields. 



It is realized for jV'*^'^''! ^ IV'^^'l in strong magnetic 



(1) 

Ae =2n. 




FIG. 17: (Color online). Plasma of fractional vortices in the 
N — 2 model in the regime |V''^''| ^ IV'*'^'! high tempera- 
tures. 




TEMPERATURE OF 
"DECOMPOSITION" OF 
COMPOSITE VORTICES 



FLUX LINE LAHICE 
MELTING 



FIG. 18: (Color online). A schematic phase diagram of dif- 
ferent phases of vortex matter and phase transition lines in 
the N = 2 model in the regime 7^ |^/''^'|. At tempera- 

tures Tm, Tm, and Tslm the melting of the composite vortex 
lattice, the sublattice of heavy vortices and the sublattice of 
the light vortices occurs, respectively. At Tlp the compos- 
ite vortices decompose. The temperatures Tl and TfJ denote 
temperatures where a phase transition via a proliferation of 
vortex loops would take place in the absence of a magnetic 
field in models with bare phase stiffnesses |7/)'"'|(T) equal to 
IV''"' |(S, T) (where a = 1, 2), if the effect of a magnetic were 
to be taken into account only via the depletion of the modulus 
of the order parameter. 



E. Vortex line plasma in the N = 2 model 



F. Physical interpretation of the external 
field-induced phases of the N = 2 model 



If the temperature is raised either at strong or 
weak magnetic fields, a situation arises where all field- 
induced composite vortices are decomposed and disor- 
dered. In addition, closed loops have proliferated2Si2S*^. 
A schematic picture of this state is shown in Fig. El 
The resulting phase diagram of the N — 2 GL model fea- 
turing the various transitions described above, is shown 
in Fig. El 



We next discuss the physical interpretation of the vari- 
ous phases that appear as a result of the above described 
vortex matter transitions. The resulting phases, which 
exhibit some quite unusual properties, come about as a 
result of the interplay between the topology of the sys- 
tem and thermal fluctuations. This is rather remarkable, 
given the three-dimensionality of the systems we con- 
sider. 
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1. Vortex lattice melting and the disappearance of 
superconductivity. 

Consider first the melting transition of an interacting 
ensemble of composite Abrikosov vortices. This phase 
transition, which is of first order—, corresponds to the 
lines T^iB) and T^iB) shown in Fig. [H It is only the 
gauge-charged mode that couples to the external field, 
while the neutral mode does not. The charged mode at 
low temperature forms an Abrikosov vortex lattice with 
a melting temperature that is suppressed with increas- 
ing magnetic field^£'^2iSfi'^. The melting temperature 
of the Abrikosov vortex lattice can be suppressed below 
the temperature where the neutral mode proliferates and 
where the composite vortex lines decompose. For A'^ = 1, 
it is known that when the Abrikosov lattice melts, super- 
conductivity is lost also along the direction of the mag- 
netic field^SiS. The situation in the N = 2 model is 
much more complex, since then there still exists a super- 
fluid mode (the gauge neutral mode) which is decoupled 
from external magnetic field. Thus, upon melting of the 
Abrikosov lattice we arrive at emergent effective neutral 
superfluidity existing in a system of charged particles^. 
This is a genuinely novel state of condensed matter, and 
moreover one which should be realizable in liquid metal- 
lic states of light atoms at in principle experimentally 
accessible pressures in the range of 400GPa!22i. 

So in the absence of an external magnetic field, the 
system thermally excites only fractional flux vortices in 
the forms of loops, with phase windings only in individ- 
ual condensates, and these fluctuations are responsible 
for critical properties. In contrast, the purely charged 
vortices (i.e. the composite one-flux-quantum vortices 
with no neutral super-flow) are not relevant in the ab- 
sence of external field and the system is either a super- 
conductor (below Tc2) or a superconductor with neutral 
mode (below Td). Thus, the effect of a sufficiently 
strong magnetic field essentially inverts the temperatures 
of the transitions by melting the lattice of charged modes 
at r^(_B) while leaving neutral modes intact. 

The phase transition from a superconducting super- 
fluid phase where the neutral mode is superfluid and 
the Abrikosov vortex lattice is intact such that longitu- 
dinal superconductivity (parallel to the magnetic field) 
exists^SiS, to a metallic superfluid phase where the sys- 
tem is superfiuid, but longitudinal superconductivity is 
lost due to the melting of the vortex lattice, can be 
mapped onto a lattice melting transition in the = 1 
model , because it is governed only by composite vor- 
tices and neutral modes are not involved. Thus it is a 
first order phase transition^Si^. 



2. Decomposition. The disappearance of superfluidity. 

Analogously, the physical meaning of the sublattice 
melting transition Tslm (see Fig. I18|) is a transition 
from a superconducting superfluid to an ordinary one- 



gap superconductor, because a disordering of the phase 
O^-^^f destroys the massless neutral boson associated with 
the gauge invariant phase difference 0*^^^ — 9'-^\ 

If we heat the system further, the ordinary supercon- 
ductivity will disappear via disordering of the phase 
when we reach the melting transition of the remaining 
sublattice of "heavy" vortices at T^. 

The system features one more phase transition. That 
is a transition from the metallic superfluid to a normal 
fluid, which has a purely topological origin. That is, from 
the vortex matter point of view, this manifests itself as 
a decomposition of a liquid of composite vortices to a 
"plasma" of individual vortices at the characteristic tem- 
perature Tlp, and such a transition has no counterpart 
in an A'^ = 1 superconductor. A schematic diagram of 
the resulting physical phases is shown in Fig. 1191 



3. A direct SSF NSF transition 

We note also the possibility of an existence of a phase 
transition directly from a superconducting superfluid 
(SSF) to a metallic normal fluid (NF), shown in Fig. 
1201 What is remarkable is that this resembles (while 
indeed being a different type of transition) the type of di- 
rect phase transition from a low temperature phase with 
Higgs mass and superfiuid density of the neutral mode, to 
a phase with zero Higgs mass and zero helicity modulus 
of the neutral mode that we would find in zero magnetic 
field when the bare phase stiffnesses of Eq. ^ are equal, 
i.e. IV'^^^I = l^/'^^^l The SSF phase features one 

massless dual Higgs photon, and one massive Higgs pho- 
ton while the NF phase features one massless photon and 
one massive dual Higgs photon. These phases are there- 
fore self-dual, in analogy to the situation encountered 
when Eq. JQ) is viewed as a quantum anti-ferromagnet 
with easy-plane anisotropy^. However, there is a sig- 
nificant difference. The novel critical point encountered 
in the case of the quantum antiferromagnet was a result 
of a superposition of a 3Dxy and an inverted SDxy crit- 
ical point. In the 2-fiavor London model in finite mag- 
netic field, the crossing point is a superposition between 
a SDxy critical line and first order phase transition line, 
so it should have a different character than the self-dual 
critical point discussed in Ref. 0. It can, however, not 
be a 3T)xy or inverted 3Dxy critical point, since neither 
of these are critical points between self-dual phases. 

Thus, both for high and low magnetic fields we have 
genuinely novel physics in that for a three dimensional 
system, i) a critical phenomenon takes place inside ei- 
ther the vortex liquid phase (high magnetic fields) or the 
vortex lattice phase (low magnetic fields), and ii) we have 
three new equilibrium states which have no counterpart 
in the A^ = 1 case. Moreover, by the discussion given 
above, it is clear that the crossing point between the vor- 
tex lattice melting line and the neutral mode prolifera- 
tion line warrants further study. This is best left for a 
separate computational analysis. 
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FIG. 19: (Color online). A schematic phase diagram of 
physical states appearing in the N = 2 model in the regime 
It/j'^'I 7^ IV'^'^'l a consequence of vortex matter phase tran- 
sitions. Increasing the magnetic field suppresses the melt- 
ing transition of the composite vortex lattice formed by the 
charged mode below the proliferation line for the neutral 
mode, which does not couple to magnetic field. In the ab- 
sence of disorder (pinning of vortices), superconductivity only 
remains along the direction of the magnetic field, provided 
that the vortex system remains in a lattice phase. When 
the composite vortex lattice melts, the system looses ability 
to carry dissipation-less charge currents, but at large enough 
magnetic field, the neutral mode should still be superfluid 
above the melting temperatura22i. Thus, we have a first or- 
der phase transition from a superconducting superfluid to a 
metallic superfluid. The neutral mode proliferates through a 
second order phase transition in the 3Dxy universality class. 
Therefore, at large enough magnetic fields, a 3Dxy anomaly 
in the specific heat should appear inside the vortex liquid 
phase. The separation between the first order specific heat 
anomaly due to vortex lattice melting and the 3T>xy anomaly 
due to loop-proliferation should increase with increasing mag- 
netic field. At low magnetic fields one has another phase 
transition inside the Abrikosov vortex lattice phase, from a 
superconducting superfluid to a one-gap ("ordinary") super- 
conducting state. 



G. Monte-Carlo results, finite magnetic field, N = 2 

We now present large-scale Monte Carlo results for the 
case TV = 2 in finite magnetic field at low magnetic 
fields when the temperature is varied^i. We consider 
the model based on Eq. for = 2 on an L'^ lat- 
tice (with L up 96) with periodic boundary conditions 
for couphng constants = 0.2, — 2, and 

e2 = 1/10. The ratio j-^^^^ | VIV"^^^ P = 10 brings out one 
second order phase transition at Tsi,m{B) in the SDxy 
universality class well below the melting temperature 



of the vortex lattice. In LMH |2/|^(i) | 



2000, 



A direct SSF to NF phase transition 
of a nonstandard universality class 




\1G-SC 



FIG. 20: (Color online). A direct phase transition from SSF 
to NF phase (shown with red arrow). 



but the physical picture remains. For real estimates of 
TsLM and in LMH, see Ref. H The Metropolis algo- 
rithm with local updating is used in combination with 
Ferrenberg-Swendsen reweighting. The external mag- 
netic field B studied is = Bv = 0, = 27r/32, 
thus there are 32 plaquettes in the (x, ?/)-plane per flux- 
quantum. This is imposed by splitting the gauge field 
into a static part Aq and a fiuctuating part Afluct- The 
former is kept fixed to (Ag, A^(r), Ag) = (0,27ra;/,0) 
where / = 1/32 is the magnetic filling fraction, on top of 
which the latter field is free to fluctuate. Together with 
periodic boundary conditions on Afiuct, the constraint 
^p(Ao -f Afluct = 27r/L^, where C is a contour en- 
closing the system in the (x,y)-plane, is ensured. It is 
imperative to fluctuate A, otherwise type-1 and type-2 
vortices do not interacliiSiiS. To investigate the transi- 
tion at TsLM we have performed finite size scaling (FSS) 
of the third moment of the action. The simulations are 
done by using vortices directly-'^^ , but with a finite mag- 
netic induction B^ — 27r/32. 

We compute the specific heat Cv and the third mo- 
ment of the action. To probe the structural order 
of the vortex system we compute the planar struc- 
ture function S''^"'(k^) of the local vorticity n("^(r) — 
(V X [V6i(°) - e A] ) /2n, given by 



1 



(/L3)2 



IE 



(r) e 



ikj_-rj_ |2\ 



(57) 



where r runs over dual lattice sites and is perpendic- 
ular to B. This function will exhibit sharp peaks for 
the characteristic Bragg vectors K of the type-a vor- 
tex lattice and will feature a ring-structure in its cor- 
responding liquid of type-a vortices. The signature of 
vortex sub-lattice melting will be a transition from a 
six-fold symmetric Bragg-peak structure to a ring struc- 
ture in 5^^^(K) while the peak structure remains intact 
in S^^\K.). Furthermore, we compute the vortex co- 
centricity Nco of type-1 and type-2 vortices, defined as 

iVco 



Nco - Nco, where 



,(2) 



nL"'(r),±ni^'(r) 



(r)| 



(58) 



where Sij 
sidering iV, 



is the Kronecker-delta. The reason for con- 
■n is that we then eliminate the effect of ran- 
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dom overlap of vortices in the high-temperature phase 
T > TgLM due to vortex-loop proliferation, and focus on 
the compositeness of field-induced vortices. 

The quantity N^o is the fraction of type- 2 vortex seg- 
ments that are co-centered with type-1 vortices, provid- 
ing a measure of the extent to which vortices of type-1 
and type-2 form a composite vortex system. Hence, it 
probes the splitting processes visualized in Fig. The 
results are shown in Fig. |^ 
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FIG. 21: (Color online). MC results for = 2 jV'^^^P = 0.2, 
IV'^'I^ = 2, and e = 1/VW- Panel a: Cv (black) and iVco 
(green). The Cv-anomaly at Tslm ~ 0.37, where type-1 vor- 
tices proliferate, matches the point at which A^'co drops to 
zero. Thus, type-1 vortices are torn off type-2 vortices. The 
remnant of the zero-field anomaly in Cy is seen as a hump at 
T ~ 3.6. Panel b: 5<^'(K) (red) and ^(^'(K) (blue) for the 
particular Bragg vector K — (tt/I, — 7r/4). 5''-"'^-'(K) vanishes 
continuously at Tslm, while ^'^'(K) vanishes discontinuously 
at = 2.34. Panel c: FSS plots of the M3 from which 
the exponents a = —0.02 ± 0.05 and 1/ — 0.67 ± 0.03 is ex- 
tracted, showing that the sub-lattice melting is a SDxy phase 
transition. Panels d, e, f, and g: plots of S'^^-'(k^) for the 
temperatures Td = 0.35, Te = 0.4, Tf = 1.66, and Tg = 2.85, 
respectively. At Td, Te, and Tf, the vortex lattice remains 
intact. The vortex lattice melts at to give a vortex liquid 
ring pattern at Tg. 

At Tslm, Cy has a pronounced peak associated with 
the SDxy transition, and a broader less pronounced peak 
which is the finite field remnant of the zero-field in- 
verted 3Dxy transition^. Scaling of M3 at Tslm shown 
in the inset c in Fig. |^ yields the critical exponents 



a = -0.02 ± 0.05 and v = 0.67 ± 0.03 in agreement 
with the SDxy universality class. A novel result is that 
S'*^^-'(K) vanishes continuously as the temperature ap- 
proaches Tslm from below, precisely the hallmark of the 
decomposition transition that separates the two types of 
vortex states depicted in Figs. El and El A related 
feature is the vanishing of iVco at Tslm as a function of 
temperature, discussed in detail below. The first-order 
melting transition takes place at Tj^, where S^'^^{'K) van- 
ishes discontinuously. This is the temperature at which 
the translational invariance is restored through melting 
of the type-2 vortex lattice. In the temperature inter- 
val T < Tslm the system features superconductivity 
and superfiuidity simultaneousljiSS, since there is long- 
range order both in the charged and the neutral vortex 
modes. In the temperature interval Tslm < T < T^ 
long-range order in the neutral mode is destroyed by 
loop-proliferation of type-1 vortices, hence superfluidity 
is lost^'^. However, longitudinal one-component super- 
conductivity is retained along the direction of the exter- 
nal magnetic field. For T > T^ superconductivity is also 
lost, hence this is the normal metallic state, which is a 
two-component vortex liquid. 

The most unusual and surprising feature is the con- 
tinuous variation of S'(^^(K) with temperature, even at 
Tslm where it vanishes. The explanation for this is the 
proliferation of type-1 vortices (which destroys the neu- 
tral superfluid mode) in the background of a composite 
vortex lattice, which the type-1 vortices essentially do 
not see, cf. Fig. {22 As far as the composite neutral 
Bose field 0^^^ — 6''^^^ is concerned, it is precisely as if the 
composite vortex lattice were not present at all. Hence, 
S'*^^-'(K) vanishes for a completely different reason than 
S'*^^)(K), namely due to critical fluctuations, i.e. vortex- 
loop proliferation in the condensate component with low- 
est bare stiffness. Such a phase transition does not com- 
pletely restore broken translational invariance associated 
with a vortex lattice, since for the type-2 vortices quite 
remarkably, the vortex lattice order survives the decom- 
position transition, due to interaction between heavy vor- 
tices mediated by charged modes. The vanishing of Nco 
is particularly interesting, and finds a natural explana- 
tion within the framework of the above discussion. That 
is, for T <C Tslm, we have A^co ~ 1, so the vortex system 
consists practically exclusively of composite vortices. As 
the temperature increases, thermal fluctuations induce 
excursions such as those illustrated in in Fig. which 
reduces from its low-temperature value, reaching a 
minimum at Tslm and then increase for T > Tslm- 

We may view the splitting process as a type-1 closed 
vortex loop superposed on a vortex lattice of (slightly) 
fluctuating composite vortices. An important point to 
notice is that a type-a vortex does not interact with a 
composite vortex by means of a neutral mode. This fol- 
lows from a topological argument that two split branches 
will feature nontrivial winding in the composite neu- 
tral field 6'^^^ — 6*'^', while a composite vortex line does 
not. Hence, the splitting transition may be viewed as a 
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type-1 vortex loop-proliferation in a neutral superfluid. 
This is illustrated in Fig. Thus, we may utilize 
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FIG. 22: (Color online). Detailed illustration of the low- 
temperature thermal fluctuations in a vortex lattice of com- 
posite vortices. A local excursion of a type-1 vortex away 
from the composite vortex lattice may be viewed as a type-1 
bound vortex loop superposed on the composite vortex lat- 
tice. The composite vortex line does not interact with a vortex 
with nontrivial winding in A7 = A(0^^' — 0*^^'). A splitting 
of the composite vortex lattice may be thus viewed as a zero- 
/je/rf vortex-loop proliferation of type-1 vortices; a 'iDxy phase 
transition universality— 

the well-known results for the critical properties of the 
SDxy model for neutral superfluids described as a vortex- 
loop proliferation^Si^Siii. This "vortex sublattice melt- 
ing" phase transition is therefore in the 3Dxy universal- 
ity class^Si^SiiZ, not a first order melting transition. The 
resulting phase is one where superfluidity is lost and lon- 
gitudinal superconductivity retained in the component 
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Conversely, remains essentially zero until Tslm, 
thereafter increasing monotonically. For temperatures 
above, but close to Tslm, fluctuations in vortices orig- 
inating in AS*^^-' are still small, so the variations in 
-^co — ~ ^co reflect thermal fluctuations in vortices 
originating in A6'^^^. The increase of means that 
type-1 vortex loops are thermally generated, and thus 
tend to randomly overlap more with the moderately fluc- 
tuating type-2 vortices. At their first order melting tran- 
sition, type-2 vortices fluctuate only slightly. Thus, the 
vanishing of N^o above Tslm reflects the increase in the 
density of thermally generated type-1 vortex loops in the 
background of a slightly fluctuating type-2 vortex lattice. 



H. Graphical representation of phase disordering 
transitions in the N = 2 model. 



In Fig. 1221 we present a schematic picture of configura- 
tion of the order parameters phases 9^^^ and 9^"^^ in var- 
ious points in physical space, when vortex matter drives 
the system into one of the above discussed superconduct- 
ing and superfluid states. 
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FIG. 23: (Color online). Phases of the order parameters in 
the various states for A'^ = 2. In the upper left panel, both S'^' 
and 61(2' are ordered, this is the superconducting superfluid 
state. In the upper right panel, neither of the phases S'^' and 
S'^' are ordered, however the combination S'^' — S'^' exhibits 
long-range order, this is the metallic superfluid state. In the 
lower left panel, S^^' is disordered and 6'^' is ordered. In 
this case, the neutral superfluid mode is destroyed and we 
are left with one charged superconducting mode, this is the 
analog of the one-gap superconducting state. In the lower 
right panel, neither of the phases 6'^' and ^'2' are ordered and 
the combination — S'^' does not exhibit long-range order, 
this is the metallic normal fluid state. The states illustrated 
in the upper left and lower right and left panels exist at zero 
as well as finite magnetic fields. The state illustrated in the 
upper right panel only exists at finite magnetic fields. 



an > 2 system is subjected to an external magnetic 
field. The new features are due to the fact that we have 
more than one neutral vortex mode, and that a vortex 
with phase winding in any single phase field 9^°^^ will ex- 
cite A^ — 1 neutral modes (see Appendix IXjl . We consider 
first the case A^ = 3, followed by the case A^ = 4. En 
route we introduce the useful concept of "color charge" 
which facilitates a discussion of the universality class of 
the phase transitions that occur in multi-component su- 
perconductors in external magnetic flcld when composite 
vortices decompose due to thermal fluctuations. 



A. Decomposition transitions, A^ = 3 



VIII. THE A > 2 MODEL IN EXTERNAL 
MAGNETIC FIELD 

We next consider the novel features that are encoun- 
tered, compared to the A' = 1- and N = 2 cases, when 



We stress that the system is not mapped onto a 
U{l)xU{l) neutral system because the neutral modes re- 
main topologically coupled, as a consequence of multiple 
connectedness of space introduced by the vortex core. It 
manifests itself in the fact that any single phase variable 
9''"^; a g [1, 2, 3] excites two neutral modes, as illustrated 
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in detail below. 

We introduce bare phase stiffnesses for the neutral 
modes in Eq. (|53|l as follows 



|^(1)|2|^(2)|2 
^2 ' 

|^(1)|2|^(3)|2 
^2 ■ 

|^(2)|2|^(3)|2 
1(2 ■ 



(59) 



Hence, a vortex with phase windings (Af?*^^-* = 
2tt,A9^'^^ = 0,A9''^^ = 0), can be mapped onto two co- 
centered vortices in a two-component neutral superfluid 
with bare stiffnesses J^^ and J^^. Thus, at a distance 
larger than the penetration length, such a vortex inter- 
acts with a vortex (A6'(i) = -27r, A6i(2) = 0, A6'(3) = 0) 
like two vortices in a neutral superfluid with bare phase 
stiffness J = J^^ + J^^. 

Intra-vortex interaction, e.g. of the vortex (A^'^' = 
27r,A6'(2) = 0,A6'(3) = 0) with a vortex {A9^^^ = 
0, A6'(2) = 27r, A0(3) ^ q) or with a vortex {A9^^^ = 
0, Ae'f^) = 0, A6'('^' = 27r) is more comphcated. It can 
most conveniently be described by introduction of the 
"color charge" concept, which we explain in the next sub- 
section, IVIII Bl 

First, however, we observe that only a composite vor- 
tex (A6'(i) = 27r, A6l(2) = 2tt, A6'(3) = 27r) has finite en- 
ergy. The key feature of a system with \ ^ 
IV'^'^-'I, is that the three elementary constituent vortices 
are bound with different strength to such a composite 
vortex. For example, when < IV'^^^I ^ I'^^^^L the 

neutral modes excited by a vortex (A6'*^^^ — 2tt, A6'^^^ — 
O,A0(''') = 0) have bare phase stiffnesses (J^^, J") < 
J^3. This in turn implies that in the composite vortex 
(A6l(i) = 27r,A6'(2) = 27r, A6i(3) = 27r), the constituent 
elementary vortex (A6'(i) = 27r, A^^^) = 0, A6l(3) = 0) 
is most loosely bound. Thus, in contrast to the N — 2 
case, the effect of thermal fluctuations for = 3 is a 
two step transition. In the general TV case the process of 
stripping a composite vortex of its N constituent vortices 
is an TV — 1-step process occurring successively, starting 
at Tci and progressing up through Tc2 up to TcAf-i, at 
which point the vortex system is fully decomposed. 

For = 3, at a low temperature determined by the 
smallest bare phase stiffness li'^^^l and by J^^ and J^^, 
there should therefore take place a partial decomposi- 
tion of the vortex (A6l(i) = 27r, A9^'^^ = 2n, A0(3) ^ 27r) 
into two vortices {A9^-^^ = 2tt,A9^^^ = 0,A9^^^ ^ 0) + 
(A6l(i) = 0, A6i(2) = 27r,A6'(3) = 2tt), illustrated in Fig. 
131 Then, upon increasing the temperature there should 
take place a phase transition, also illustrated in Fig. 1241 
into a fully decomposed state deflned by the phase wind- 
ings (A6l(i) = 27r,A6i(2) = 0, A6'(3) = 0) -f (A6'(i) = 
0,A6'(2) = 27r,A6i(3) = 2t:) -> (A6l(i) = 27r, A6'(2) = 
0,A6i(3) = 0) -I- (A6l(i) = 0,A6'(2) = 27r,A6i(3) = 0) + 
(A0(i) = 0, A6l(2) = 0, A6'(3) = 27r). 

We also stress that apart from the neutral mode 
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FIG. 24: (Color online). A two-step decomposition tran- 
sition in the A*' = 3 model in external magnetic field. The 
"lightest" vortex component, originating in the order param- 
eter component with the lowest bare phase stiffness, tears 
itself loose from the composite Abrikosov vortex of the two 
stiffer order parameter components at Tspi. At a higher tem- 
perature Tlpi, the "next-to-lightest" vortex component (green 
vortex), originating in the order parameter component with 
the next-to-lowest bare phase stiffness, tears itself loose from 
the vortex of the stiffest order parameter component in the 
background of proliferated vortices originating in non-trivial 
phase-windings of the phase with lowest phase stiffness (red 
vortex). 



(6)(2) _ g)(3))^ .^^i^ich provides an attractive interaction 
for the vortex pair {A9^^^ ^ 0, A6l(2) = 27r, A6l(3) = 
0) + (A6i(i) = 0,A6i(2) = 0,A6l(3) = 27r), the flrst 
of these vortices excites the neutral mode (6'^^^ — 9'^^^) 
which consists of oppositely directed currents in the con- 
densates "i'^'^ and ^'o^''. The second vortex excites the 
mode {9^^'> — 9^^'>) which is associated with oppositely 



directed currents in condensates '^'^^ and ^'q'^''. These 
modes, apart from giving the pair (A6'^^^ — 0, A6'^^^ = 
27r, A6'(3) =0) + (A6'(i) = 0, A6'(2) = 0, A6'(3) = 27r) log- 
arithmically divergent energy per unit lengt h^^i^^ , also 
yield some repulsive interaction in this pair, because both 
the modes (6''^' — 6*'^') and {9^^^— 9^^^) feature unscreened 
currents of condensate ^tp^^. However, such a repulsive 
interaction in this pair is negligibly small in the consid- 
ered regime IV'^^^I ^ IV'^^''! ^ IV'^^''Ij compared to the 
interactions mediated by the mode {9^^^ — 9^^^). 



,(3) 



B. Color electric charge 

Formally, the partial decomposition process can be 
described by introducing the concept of "color electric 
charges". That is, we may introduce e.g. "-|-red", 
"-|-green" and "+blue" charges associated with 27r wind- 
ings in {9^^^ - 6'(2)), (6»(i) - 0(3)) and (6l(2) - 0(3)), re- 
spectively (see Tab. 1^). If we have a — 27r winding in 
(6/(1) _ 0(2)), (0(1) - 0(3)) or (0(2) - 0(3)), that would cor- 
respond to " —red" , "—green" and "—blue" color electric 
charges, respectively. We stress once more that in order 
to preserve single- valuedness of the order parameters, the 
±27r gains in phase differences may only come as ±27r 
gains in individual phases. For example if we would have 
(A0(i) = 37r/4, A0(2) = -57r/4) then (0(i) - 0(2)) would 
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TABLE 11: A color charge is defined as a ±2tt winding in the 
mode (6l(°' - e*'"') with the following mapping for N = 3. 



"Vortex"-"Color Charge" duality 



Mode: (e*^' - f'^)) 


(0(1) _ 0(3)) 


(0(2) _ 0(3)) 


Color charge: red 


green 


blue 



change by 2it. However, such a configuration would be 
unphysical because individual order parameters ^'q^^ and 

(2) 

^'q would loose their single- valuedness. Then, a vortex 
{Ae^^^ = 27r, A6l(2) = 0, AO'-^^ = 0) which excites two 
neutral modes associated with and (6''-^^— 6''^') 

with stiffnesses J^^ and J^^, respectively, may be viewed 
as a color charged string with color charge "-|-red" and 
"-l-green" . 

The regime \^^^^\ < IV'^^^I < IV'^^^L i-e- when 
J^^ <C J^^ <C J^^, corresponds to the situation where 
red electric charges are much weaker than green charges, 
which in turn are much weaker than the blue charges. 
The blue charge then dominates the binding of the vor- 
tices (A6i(i) = 0,A6'(2) = 27r,A6l(3) = 0) and (AO^^^ = 
0, A6I(2) = 0, A6'(3) = 27r). The tightly bound compos- 
ite vortex {Ae'^^^ = 0, A6i(2) = 27r,A6l(3) = 27r) then 
has electric charge (—red, —green) which loosely binds 
it with (-|-red, -|-green) color charged vortex (A6'^^^ — 
27r, A6l(^^ = 0,A6'(3' = 0) into a color charge neu- 
tral finite energy one-fiux-quantum vortex (A0^^^ = 
27r, A6'(2) = 27r,A6l(3) = 27r). 

In Fig. we illustrate how to connect the vortex 
picture to the picture of color charges, for the case N = 3. 
For N = 3, each type of vortex is a bound state of two 
color charges. 

A schematic picture of the low-temperature composite 
vortex lattice phase, the partial decomposition transition 
in the color electric charge representation, and the com- 
plete decomposition transition, are given in Figs. 1261 1281 
and El 

For the purposes of determining the universality class 
of this partial decomposition, which is a genuine phase 
transition, it may also be viewed as follows. The com- 
pletely intact composite vortex is color charge neutral in 
the sense that it contains a positive and negative elec- 
tric charge of each color. The fluctuation snapshot in 
the left part of Fig. [31 can be viewed as a completely 
intact composite vortex line with a small type-1 vortex 
loop superimposed on it as shown on Fig. |23 The 
type-1 vortex loop which carries a red and green elec- 
tric charges does not interact with the composite vor- 
tex which is color charge neutral. This can also be 
seen by examining the vortex interaction matrix given 
in Eq. H12|l . That is, when we add up the neutral- 
mode- mediated interactions between the type-1 vortex 
with all three type-1, 2, 3 vortices in the composite vor- 
tex (A6l(i) = 27r,A6'(2) = 2n,Ae'^^^ = 2tt) they add to 
zero. The situation in the right part of Fig. [21 is topo- 
logically equivalent to a completely intact composite vor- 
tex line superimposed with one segment of an unbound 
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FIG. 25: (Color online). The connection between vortex 
illustrations and bound states of color charges, for the case 
= 3. Each type of vortex is a bound state of 2 color charges. 
The radius of each black circle is for graphical convenience is 
used to differentiate between "heavy" and "light" vortices. 
The order parameter component with the lowest bare phase 
stiffness is taken to be the vortex with the "smallest diame- 
ter". A vortex originating in a non-trivial phase- winding in 
is denoted a type-a vortex. The color of the vortex on 
the top of the Figure should not be confused with the color of 
the electric charges in the "dual" charge representation. 



vortex loop. Therefore the transition may be viewed as 
an Onsager vortex loop proliferation transitiort^^'^^^ of 
type-1 vortex loops in the background of a color charge 
neutral composite vortex latticsmSiM^. That is because 
type-1 vortex loops, from the point of view of superfluid 
modes, cannot "see" color charge neutral vortex lines, 
so this is equivalent to a type-1 vortex- loop proliferation 
transition in the complete absence of a composite color 
charge neutral vortex lattice. Since these vortices excite 
neutral modes, this transition, which is the first stage 
in decomposing a color charge neutral composite vortex 
line, is a vortex loop proliferation transition in the 3T)xy 
universality clas322i2i. In the color charge representation 
given in Fig. 1281 this also means that the first-stage par- 
tial decomposition transition of three charged fluctuating 
line objects with different charges (but such that the alge- 
braic sum of their charges add up to zero) is also a 2-color 
metal-insulator phase transition in the 3T>xy universality 
class, involving flexible color line charges. 

The usefulness of the color charge representation be- 
comes particularly clear when we go on to describe the 
second stage of the decomposition transition, illustrated 
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FIG. 26: (Color online). A color charge representation for 
= 3 of the composite vortex lattice phase illustrated also 
in the left part of Fig. 1241 This low-temperature phase, 
where the fluctuations of the composite vortex lines only in- 
volve small excursions of each of the constituent vortex lines 
away from the main composite object, may be viewed as a 
3-color dielectric phase. The main composite vortex line is 
"anchored" on the thickest vortex. We assume we are at tem- 
peratures and fields such that the thickest vortex line does 
not undergo a melting transition. Each constituent vortex 
may be viewed as a bound state of certain combinations of 
color charges. The composite vortex line may, on the other 
hand, be viewed as bound states of ±red, ±green, and ±blue 
charges, as indicated by the dotted ellipses. This is a three- 
color dielectric "insulating" phase. We strongly emphasize 
that the above illustration is meant to illustrate what the sit- 
uation is in a typical cross section along the lines, which are 
not rigid straight vortex lines. 




FIG. 27: (Color online). A local thermally driven detach- 
ment of the type-1 vortex line (red color) from the composite 
vortex line can be viewed as a superposition of thermally cre- 
ated closed type-1 vortex loop on a completely composite vor- 
tex line. Both processes are topologically equivalent because 
the vorticity of the type-1 constituent vortex in the compos- 
ite vortex line is exactly canceled by a superimposed type-1 
vortex ring with the opposite vorticity. We stress that if left 
segment of a vortex loop has a counterclockwise vorticity then 
the opposite right segment has a clockwise vorticity. 



in Fig. 1241 This transition may be viewed as a prolifera- 
tion of type-2 vortex loops in the background of hberated 
type-1 vortex loops, all superimposed on a background 
composite vortex lattice. The situation therefore is more 
complicated than in the first stage illustrated in Fig. 1241 
since that was a proliferation of type-1 loops in vacuum. 

However, let us view this transition in the color charge 
picture, illustrated in going from Fig. [2H1 to Fig. 1^ 
This is a metal-insulator transition for the blue-charge 




FIG. 28: (Color online). A color charge representation for 
A = 3 of the partial decomposition transition In the color 
charge representation, this may be viewed as a red-green 
dielectric-metal transition, while the blue dielectric phase re- 
mains intact. As explained in the text, this partial 2-color 
metal-insulator transition involving fluctuating line charges, 
is in the SDxy universality class. The dotted ellipse indicates 
which color is involved in forming the remaining dielectric 
phase. 



sector in the background of coexistent red and green 
metallic phases. However, red and green charges can- 
not screen blue charges, while these charges eliminate 
the neutral modes associated with them (that is the ones 
with bare stiffnesses J^^ and J^^). Therefore, this is 
a metal-insulator transition for the sector of the blue 
charges, while red and green charges screen themselves 
and do not affect this transition. 

In the color charge representation given in Fig. this 
also means that the second-stage decomposition transi- 
tion depicted in Fig. |^is a 1-color metal- insulator tran- 
sition in the 3Dxy universality class, involving flexible 
color charged strings. 



C. Decomposition transitions, A = 4 

In Figs. 1201 El OH and |2S1 we illustrate the partial 
decomposition transitions for the case = 4 in the color 
charge picture. That is, we introduce as in the case A^ = 3 
case "± red" , "± green" and "± blue" charges associated 
with ±27r windings in (6'(i) -6l(2)), (6»(i) -6»(3)) and (6'(2) - 
f?*^^^), respectively. In addition, we introduce "± yellow", 
"± violet" , and "± orange" charges associated with ±27r 
windings in (^(i) - 6l(4)), (6»(2) _ gi^) ) and (6'(3) - 6»(4)), 
respectively (see Tab. IIII|I . Therefore, the case A^ = 4 
features one new aspect which was absent in the case 
A^ = 3, namely that for A^ = 4 we need more color charges 
than number of order parameter components in order 
to completely cover all the possible ways that a neutral 
mode (6'(") -6'('')) can be excited. Each type-a vortex is 
a bound state of three color charges, as indicated in Fig. 

EOl 
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FIG. 29: (Color online). A color charge representation of the 
complete decomposition transition also illustrated in Fig. 1241 
In the color charge representation, this may be viewed as a 
blue dielectric-metal transition, in the background of a red- 
green metallic phase. As explained in the text, this complete 
"1-color" metal-insulator transition involving fluctuating line 
charges is in the same universality class as the partial decom- 
position transition, namely the 3T)xy universality class. 

The low-temperature phase is a 4-composite color 
charge neutral vortex system, which may alternatively 
be viewed as a 6-color dielectric phase, Fig. QUI 
The first-stage partial decomposition involves a type- 
1 vortex tearing itself off tlie 2-composite vortex, i.e. 
a vortex loop proliferation of type- 1 loops carrying 
-|-red, -l-green, -I- yellow color charges in the background 
of color charge neutral objects. So this is a phase tran- 
sition in the SDxy universality class, ft may alterna- 
tively be viewed as a 3-color metal insulator transition in 
the red, green, yellow line-charge sectors. Fig. |^ leav- 
ing a 3-color (violet, blue, orange) dielectric phase. The 
second-stage decomposition process involves a type-2 vor- 
tex tearing itself off a 3-composite vortex in the bacfc- 
ground of a system of proliferated type-1 loops. Due to 
screening of red, green and yellow charges, this transition 
may be viewed in a simplified manner, ft may be consid- 
ered as the first-stage decomposition in a iV = 3 system 
consisting of type-2, type-3, and type-4 vortices, involv- 
ing violet, blue, and orange charges. This we have already 
argued is a SDxy (type-2) vortex loop proliferation tran- 
sition, when we considered the A*" = 3 case. Alternatively, 
it may be viewed as a 2-color metal-insulator transition 
in the violet and blue line-charge sectors. Fig. leav- 
ing a I-color (orange) dielectric phase. The third-stage 
decomposition may be viewed as a type-3 vortex loop 
proliferation in the background of type-1 and type-2 pro- 
liferated vortex lines. Due to screening of violet and blue 
charges this may be viewed in a simplified manner, ft 
may be considered as a vortex loop proliferation of loops 
carrying orange charges in the background of an orange- 
neutral vortex lattice, or vacuum. This is a vortex loop 



TABLE III: A color charge is defined as a ±2tt winding in 
the mode {6^"^ — S'''') with the following mapping for A'^ = 4. 

Mode: (e(i) -e(2)) ^^wZ^) ^^wT^ 

Color charge: red green blue 

m^: (0W -0W) (e(2)_g(4)) (0(3) _eW) 
Color charge: yellow violet orange 



proliferation in the 3Dxy universality class^LSi. Alter- 
natively, it may be viewed as a 1-color (orange) metal- 
insulator transition, (see Fig. I33|l . leaving the 6-color 
dielectric phase completely destroyed. 
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FIG. 30: (Color online). A color charge representation for 
A'^ = 4 of the composite vortex lattice phase. This low- 
temperature phase, where the fluctuations of the composite 
vortex lines only involve small excursions of each of the con- 
stituent vortex lines away from the main composite object, 
may be viewed as a "6-color dielectric" phase. Moreover, for 
A = 4, each type-a vortex in the A = 4 case is a bound state 
of 3 color charges. The dotted ellipses indicate which colors 
are involved in forming the dielectric phase. 



D. General N 

In the general A^-component case, the number of color 
charges that needs to be introduced to give an equiva- 
lent description as the above, may be counted as follows, 
starting from the third term in Eq. (0}. Each combina- 
tion 0*^"^ — 9^^"^ is given a color. We start with one phase, 
the one with the lowest bare stiffness say 9^^\ and intro- 
duce a non-trivial phase winding ±27r in this phase. This 
excites N — 1 neutral modes since A'^ — 1 gauge-invariant 
phase differences which involve can be formed, fn- 
troducing non-trivial phase-windings in the next phase 
9^^\ the one with the next-to-lowest bare stiffness say, 
will also excite A^ — 1 neutral modes, but only N — 2 new 
neutral modes. Non-trivial phase-windings in the third 
phase fl^'^-' will excite N — 3 new neutral modes, and so 
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FIG. 31: (Color online). A color charge representation for 
= 4 of the first partial decomposition transition. In the 
color charge representation, this may be viewed as a red- 
green-yellow dielectric-metal transition, while the violet-blue- 
orange dielectric phase remains intact. As explained in the 
text, this partial 3-color metal-insulator transition involving 
fluctuating line charges, is in the SDxy universality class. The 
arrows indicate which three colors are involved in the metal- 
insulator transition. The dotted ellipses indicate which colors 
are involved in forming the remaining dielectric phase. 



on. The number of different colors A^coior we will have 
to introduce for a theory with N flavors of scalar fields 
is therefore given by A^coior(^) = (A^ - 1) + (A^ - 2) + 
(TV - 3) -f • • • + 2 + 1 = NiN - l)/2, i.e.N,oioA2) = 1, 
A^coior(3) = 3, 7Veoior(4) = 6, and Arcolor(5) = 10. A 
completely composite vortex, which we denote as an iV- 
composite vortex, consists of N constituent vortices origi- 
nating in nontrivial phase windings in each of the individ- 
ual phases 6^"\ a & [1, . . . , iV]. Since a non-trivial phase 
winding in any phase 6'-^\ rj € [1, . . . , A^], excites iV — 1 
neutral modes 6'^''^ — 9'^"\ it is clear that a type-77 vortex 
may be viewed as a bound state of — 1 color charges. 
The particular combination of A^ — 1 color charges, out 
of the total collection of N{N — l)/2 color charges, that 
will enter the A^ — 1-body bound state in each vortex, will 
depend on rj. The N -composite vortex is a color charge 
neutral object. 

Small fluctuations in the A"-composite vortex may 
therefore be viewed as a dielectric insulating phase of an 
A^(A^ — l)/2-component dielectric. The first stage in the 
A^ — 1 stage decomposition process of the A^-composite 
vortex line, where a type-1 vortex tears itself off the A^- 
composite vortex, is therefore a metal-insulator transi- 
tion where A^ — 1 color charges of the A^(A^ — l)/2-colors 
dielectric system simultaneously undergo a metal insula- 
tor transition, in the 'SDxy universality class. The next 
stage, where a type-2 vortex tears itself off the remain- 
ing A^ — 1-composite vortex in the background of a sys- 
tem of proliferated type-1 loops, is phase where N — 2 
color charges simultaneously undergo a metal insulator 
transition in the 3Dxy universality class by the same ar- 
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FIG. 32: (Color online). A color charge representation for 
A = 4 of the second partial decomposition transition. In the 
color charge representation, this may be viewed as a violet- 
blue dielectric-metal transition, while the orange dielectric 
phase remains intact. As explained in the text, this partial 
2-color metal-insulator transition involving fluctuating line 
charges, is in the SDxy universality class. The arrows in- 
dicate which two colors are involved in the metal-insulator 
transition. The dotted ellipse indicates which color is involved 
in forming the remaining dielectric phase. 



gument as used for the A^ = 3 case, and so on. The 
complete decomposition of the A"-composite vortex pro- 
ceeds in Af — 1 steps of metal-insulator transitions for 
color charges, where step number 1 < Af < N — 1 may 
be viewed as either a type-A/" vortex tearing itself off 
an N — Af — 1-composite vortex line, or equivalently a 
simultaneous metal-insulator transition for N — Af new 
color charges that have not been involved in the previous 
Af — 1 metal-insulator transitions. All the A^ — 1 partial 
decomposition transitions, or metal-insulator transitions 
for color charges, are in the 3Dxy universality class. 



We should emphasize that, as follows from Eq. Q in 
the limit A^ — > 00, the strength of each of the electric 
charges goes to zero. At the same time the number of 
colors of electric charges Nc tends to infinity. From Eq. 

it follows that even in the limit N ^ 00, the energy 
binding of a type-a vortex to a color charge neutral com- 
posite vortex is finite, even though the strength of each 
individual color charge tends to zero. 
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FIG. 33: (Color online). A color charge representation for 
A'^ = 4 of the third and complete decomposition transition. 
In the color charge representation, this may be viewed as an 
orange dielectric-metal transition, while the color-dielectric 
phase is completely destroyed. As explained in the text, this 
partial 1-color metal-insulator transition involving fluctuat- 
ing line charges, is in the 3Dxy universality class. The arrow 
indicates which color is involved in the metal-insulator tran- 
sition. 



E. Graphical representation of phase disordering 
transitions in the A = 3 model. 

In Fig. 1^ we illustrate graphically the various phase 
disordering transitions and partial symmetry restorations 
discussed in the previous section. 



IX. APPLICATIONS 

In this Section, we will briefly mention some possible 
applications of the results obtained in this paper. Em- 
phasis will be on the case N — 2, but some results for the 
cases iV = 3 and even = 4 may find applications in 
mixtures of superconducting condensates in the not too 
distant future. 



A. Applications of results for N = 2 

The results for N = 2 are expected to apply to two- 
component superconductivity which could be achieved in 
metallic states of light atomsS*^, such as electronic and 
protonic condensates in liquid metallic hydrogen (LMH) 



under extreme pressure. Estimates exist for Tc2 for such 
systems, Tc2 ~ 160K^. A rough estimate for Td follows 
from the mass ratio of the electronic and protonic con- 
densate, Tci « O.IK. Hence, at Td one should observe 
an extra low-temperature SDxy specific heat anomaly, 
as well as an anomaly in the London penetration depth. 
An even more promising candidate is the the system 
CH^ where there are predictions of liquid metallic states 
at considerably lower pressures than those required to 
achieve LMH. In such multi-flavor superconductors it is 
the case IV'^^^I ^ IV'^^^I which appears most naturallyJsi. 

Here, it is appropriate to remark briefly on the micro- 
scopic origins of superconductivity in the projected liquid 
metallic phase of hydrogenSi^. The proton is four times 
lighter than a "^He atom. It is well known that '*i/e at 
normal conditions is a classic permanent liquid^ because 
of high zero-point energy and weak ordering energies. In- 
deed zero-point energies of protons in a dense environ- 
ment are also high, and at increasing compression there 
is a shift of electron density from intra-molecular regions 
to inter-molecular, and with it a progressive decline in the 
effective inter-proton attractions. Because of this there 
is also a decline of ordering energies from interactions 
relative to protonic zero-point energies. The existence 
of a melting point maximum as a function of pressure in 
hydrogen, as well as a range of densities where hydrogen 
may take up a fluid phase in its ground state was sug- 
gested in Ref. 0. Another important circumstance is that 
en route hydrogen should undergo an insulator-metal 
transition and therefore the resulting phase should be the 
liquid metallic hydrogen, a translationally invariant two- 
component fermionic liquid. There is preliminary exper- 
imental evidence that a melting point maximum may in- 
deed exist^ and it has received recent powerful backing in 
ah initio calculations^. Experimentally a 12.4 fold com- 
pression of hydrogen has already been achieved at around 
320 GPa^'^. Estimates suggest that LMH should appear 
at 13.6 fold compression at pressure in the vicinity of 400 
GPaS, whereas hydrogen alloys may exhibit metallic be- 
havior at significantly lower pressures^. A predicted key 
feature of LMH at low temperature is the coexistence of 
superconductivity of proton-proton and electron-electron 
Cooper pairs^. 

The special point = has a physical realiza- 

tion when Eq. ^ is viewed as an effective field theory of 
a quantum antiferromagnet with easy plane anisotropy, 
which facilitates a suppression of topological defects in 
the form of "hedgehog" -configurations which appear in 
0(3)-symmetric models^. (More generally, it may be 
viewed as a field theory of an 0(3) model where "hedge- 
hogs" are suppressed by some unspecified mechanism, 
not necessarily limited to easy-plane anisotropy). 

We conclude this subsection on the N — 2 case with 
a remark on how these results relate to multi-flavor elec- 
tronic condensates. To describe this case, we need to 
include a Josephson coupling between the matter field 
species. The details of how to give a vortex represen- 
tation of the iV-flavor London model in the presence 
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of inter-flavor Josephson coupling, is given in Appendix 
IE] Had a Josephson coupling term between condensate 
species been introduced in the theory for > 2, this 
would have have altered the dual theory in a completely 
non-perturbative way, and would tend to lock the phases 
of individual condensate fields to each other. (For a dual 
representation of this case, where the non-perturbative 
character of the Josephson coupling is brought out in a 
particularly clear way, see Appendix lEjl . As a result, 
the transitions we describe here would collapse to one, 
namely the charged Higgs fixed point, which is in the 
inverted 3Y)xy universality class. 



B. Applications of results for A^ = 3, 4 

Mixtures of superconducting condensates in LMH can 
be extended to include also the hydrogen-isotopes deu- 
terium and tritium.^^. Tritium is a S' = 1/2 fermion, 
so this may give rise to a superconducting condensate 
via forming spin-singlet Cooper pairs, just as in the pro- 
tonic case. Hence, our results for A^ = 3 could be ap- 
plicable to to the mixtures of liquid metallic hydrogen- 
tritium at extremely high pressures. Another possibility 
is to include deuterium as a new component. Includ- 
ing deuterium as a new component in addition to hy- 
drogen means that we have an A^ = 3-component mix- 
ture of superconducting condensates consisting of elec- 
trons, protons, and deuterons (deuterium nuclei), all of 
which in principle can undergo a metal-superconductor 
transition. Compared to the situation for N = 2, the 
situation is complicated by at least two circumstances. 
Firstly, deuterons are bosons, and secondly they have 
spin S = 1. This means that the electrons and protons 
become superconducting via forming Cooper pairs, while 
the deuterons undergo a metal-superconducting transi- 
tion via Bose-Einstein Condensation. Extending this to 
the case of having both tritium and deuterium in addi- 
tion to hydrogen, might provide a realization of the case 
A^ = 4. 



X. SUMMARY 

We have analyzed the A^-flavor London superconduc- 
tor model coupled to one gauge field with no Josephson 
coupling between the matter field components. The dual 
theory is an A^-flavor GL theory coupled to A" dual gauge 
fields, where the sum of all dual gauge fields is massive at 
all couplings. There are A^ — 1 charge-neutral superfluid 
modes and one charged superconducting mode in this 
model. We have given a prescription for how to identify 
the N — 1 neutral modes for arbitrary A^. 

For A^ = 2, a case which should apply to a supercon- 
ducting state of liquid metallic hydrogen, as well as for 
the case A^ = 3, we have performed large scale MC simu- 
lations computing i) critical exponents a and i^, ii) gauge 
field and dual gauge field correlators. Hi) the correspond- 



ing masses, and iv) critical couplings using ESS. For 
N = 2 and l^/"*-^-*! ^ |V'^^''Ij we find one low-temperature 
critical point in the 3Dxy universality class at Tci, and 
one critical point in the inverted SDxy universality class 
at Tc2 > Tel. Eor A^ = 2 with ji/'^^^l = IV'^^^I we find 
one critical point with non-SDxy values of a and v^. We 
propose these to be critical exponents in the universality 
class of a self-dual gauge theory. For A^ = 3 and all IV''-"' | 
unequal, we find two fixed points in the 3T)xy universality 
class at Tci and Tc2, and one fixed point in the inverted 
3Dxy universality class at Tc3(> Tc2 > Td). All critical 
points therefore exhibit SDxy values of a and u when 
all bare phase stiffnesses are different. In the case 

AT = 3 with IV'^^^I = \tp^^^ < IV'^^^I we find two critical 
points. The critical point at Td < T'c2 is found to be in 
the SDxy universality class, and the critical point at Tc2 
is in the inverted SDxy universality class. For the case 
A^ = 3 with equal phase stiffnesses we find one critical 
point which is non-SDxy . 

In this context, we have also noted that collapsing two 
neutral critical points in the 3Dxy universality class leads 
to a single critical point also in the SDxy universality 
class. This follows from an argument implying that col- 
lapsing any number of neutral critical points in the SDxy 
universality class leads to a single critical point also in 
the SDxy universality class. On the other hand, it ap- 
pears that collapsing N — 1 neutral critical points in the 
SDxy universality class and one charged fixed point in the 
inverted SDxy universality class leads to a single critical 
point in a universality class (which in principle could de- 
pend on A^) which is not that of the SDxy or inverted 
SDxy type. For N = 2, we may define the universality 
class as that of a 3D self-dual U{1) x U{1) gauge theory, 
while it is less clear what it is for other A^ > 3. The 
numerical values we have obtained for the critical expo- 
nents a and v the two cases i) N = 2, \^^^^ \ = IV'*'^''! and 
ii) N — S, IV''"'"-'! = IV''"^''! = IV^*"^"*! are remarkably similar, 
indicating that the values of the critical exponents are at 
most weakly dependent on A^. 

In an external magnetic field at low temperature, the 
ground state of the system is an Abrikosov vortex lat- 
tice of composite vortices. However the effect of thermal 
fluctuations alters the physics signiflcantly. We discuss 
in detail that in the low-field regime and when the bare 
stiffnesses of the condensates all differ, we find that a 
SDxy vortex sub-lattice melting transition takes place. 
Upon thermal decomposition of field induced compos- 
ite vortex lines, the constituent vortices originating in 
the condensates with lowest bare stiffness disorder, while 
the ones originating in the stiffer condensates remain ar- 
ranged in an Abrikosov vortex lattice. When such a 
transition occurs, an A^ = 2 system looses superfluid 
properties, but remains superconducting. In contrast at 
high magnetic fields, the neutral mode disappears via the 
melting transition of the lattice of composite vortices. 
This is a transition from the superconducting state to 
a non-superconducting metallic superfluid stateSS. In- 
side this novel metallic superfluid phase, at a tempera- 
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ture much lower than the zero- field metal-superconductor 
transition, we find a superfluid-normal fluid SDxy phase 
transition associated with a neutral mode-driven prolif- 
eration of vortex loops nucleating on field-induced vor- 
tex lines. These various transitions should in principle 
be detectable in flux-noise experiments. We extend the 
discussion for N > 2, where we find phase transitions 
associated with partial decomposition of composite vor- 
tices, yielding several unusual states of partial symmetry 
breakdown. The universality class and partial symme- 
try breakdowns are identified by mapping the system to 
an ensemble of electrically charged strings where for the 
A^-component system there are N{N — l)/2 replicas of 
electric charges of different color. 

The sublattice melting and partial decomposition tran- 
sitions are of purely topological origin. It involves what 
can be viewed as vortices and anti-vortices (positively 
and negatively charged objects). The existence of pos- 
itively and negatively charged objects implies that the 
physics conceptually in some sense is similar to what oc- 
curs in a Kosterlitz-Thouless (KT) transition in two di- 
mensions. In spite of being a decomposition of positively 
and negatively charged composite objects, such a transi- 
tion can be mapped onto a proliferation of vortex loops in 
a vacuum, Lea, phase transition in the SDxy universality 
class. Another principal difference between this type of 
phase transition compared to the KT transition, is that 
in the 2D KT transition, vortices and anti-vortices (posi- 
tively and negatively charged objects) are thermally gen- 



erated. This can not happen in three dimensions since 
any vortex line in 3D has an infinite energy in an infinite 
sample. In the 3D transition which we consider, the neu- 
tral bound states of the charged objects are introduced 
by an external magnetic field. Thus we deal with a novel 
type of transition which in a very unusual form involves 
concepts both from two- and three- dimensional physics. 
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APPENDIX A: IDENTIFYING CHARGED AND NEUTRAL MODES 



Here, we illustrate in detail the separation of variables in iV-component London model and the extraction of the 
composite charged and neutral modes in the London limit. The general procedure beyond the London limit for the 
case N — 2 can be found in^^. The general A-component Ginzburg-Landau model is given by the action 



f ^ 1 1 

^ = E E 2l(^ - ieA)V'(")|^ + V^({V'^")}) + -(V X A)^ 

r {a=l 



(Al) 



Here, the masses M*^"^ have been absorbed in the amplitudes for notational simplicity. Let us rewrite the first 

term in Eq. (jAl|) as follows 



N N 

-|(V - ieA)^(")|2 =. J2 i|^^"^P(V0("))2 - eA • (Ve^"))!^^^") p + ie^A^I^-^") p. 



(A2) 



The idea now is to first extract the charged mode of the system, and then identify the remaining terms as the neutral 
modes. The charged modes is the (only) linear combination of phase gradients that couples to the gauge field 

A. We first introduce the quantity 



N 



= ^ \i;^°'^\^ve^°'\ 



(A3) 



Using this, we form the combination 



N 



E iv^'" 



)|2 



V6l(") - eA 



(©-e^'^A)^ 



(A4) 
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Phases 

of the fields: 




SUPERCONDUCTING SUPERFLUID 
WITH TWO NEUTRAL AND 
ONE CHARGED MODES 



SUPERCONDUCTING SUPERFLUID 
WITH ONE NEUTRAL AND 
ONE CHARGED MODES 



METALLIC SUPERFLUID 
WITH TWO NEUTRAL AND 
NO CHARGED MODES 




FIG. 34: (Color online). A schematic plot of states in the 
A*' = 3 system. The upper left panel shows the phases of 
the condensate order parameters that are involved. The up- 
per right panel shows the state where all three phases 6'^^\ 
9^^\ and 9^^^ are ordered individually. This state is the low- 
temperature (ground) state and features one superconducting 
charged mode and two superfluid neutral modes. The mid- 
dle left panel is a phase where 6^^^ is disordered, while 6^^^ 
and e'^" are ordered. Thus, this is a state which features 
one charged superconducting mode and one neutral super- 
fluid mode. The middle right panel illustrates a state where 
aU of the phases g'^' , 6^^^ , and gt^' are individually disordered. 
However, the differences 6^^'> - O'-'^^ and f*^' - 9'-'^'> (and there- 
fore also 9^^^ —9^^^) feature long-range order. This is therefore 
a state which is normal metallic, but nevertheless features two 
neutral superfluid modes. The bottom left panel illustrates a 
state where all of the phases 9^^\ 9^'^\ and 9^-^^ are individu- 
ally disordered. Only the phase difference 9^^^ — 6*'^' exhibits 
long-range order. This is a normal metallic state featuring one 
neutral superffuid mode. The bottom right panel illustrates 
a state where all of the phases 9^^\ 9^'^\ and 9^^^ are indi- 
vidually disordered and where none of the phase differences 
0(1) - 0(2) and 6l(2) - S^^^ and 0'^) - S^^' feature long-range 
order. This is therefore a state which is normal metallic and 
normal fluid (no neutral superfluid modes). The states illus- 
trated in the upper left, middle left, and lower right panel 
exist at zero as well as finite magnetic fields. The states illus- 
trated in the middle right and lower left panels only exist at 
finite magnetic fields. 
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where ^'^ is defined in Eq. |SJ- One can check that this combination is gauge- invariant. By adding and subtracting 
@2^2vlf2^ we now express Eq. (|A2() as follows 



^ 1 1 



)|2 



a = l 

1 
2 



TV 



©2 



2^(0 "^*'A)^ 



(A5) 



The last term is identified as the charged mode. The first term on the right hand side of Eq. ljA5|l can be rewritten 
as follows 



N 



©2 



1 

1 

4^ 



N 



■-.13=1 
N 

j2 iV'^"^nv'^'^^i'(v0^"^ - v^c^))- 



a, 13=1 



(A6) 



Therefore, the action may be expressed as 



TV 

^ |^(")|2|^(/5)|2 (V(0(") 



a,/3=l 



,(A7) 



which is now a sum of one charged mode and A'^ — 1 neutral modes. The last term in Eq. (|A7(I is a sum consisting of 
N{N — l)/2 terms, involving all the color charges (0^") — 0^^)) (see Section IVlIip which are needed to account for all 
the possible ways neutral modes can be excited in the system as a consequence of multiple connectedness of physical 
space in the presence of vortices. Note how this term vanishes when = 1. 



APPENDIX B: DERIVATION OF EQS. ^ AND JT^ 

Starting from the model in the Villain approximation Eq. lO we linearize the kinetic energy terms by introducing 
N auxiliary fields v'"^ where a € 1, . . . , A with a partition function 



PA [| / Vd'^-'^ / Pv'-^) exp(-5), 
^ - E f E ^^(-^"^)' + E KA^^"^ - + 2™(")) . v(") + f (A X Af 

r \q=1 ''^'^ ' a=l 



(Bl) 



where = \^\^\'^ / M'^°'\ The Poisson summation 

formula reads 



E 



JiTzinB 



(B2) 



n= — OQ 



m— — 00 



Here, n,m e Z and B & R. We apply Eq. (|B2|) to 
Eq. ljBl|) and integrate out the integer fields n*^"' so that 
the fields v*^"' take only integer values which we denote 
v^"^. After a partial summation of X^r ^^6*'^"^ • 



>(") - 



J2t Sa=i i^^"^ A • v^"^ , where the surface terms 



are omitted, we may integrate out the phase fields 
This integration produces the local constraints 



A- 



>(«) _ 



= 0. 



(B3) 



To fulfill this constraint we let v*^"-' = A x h^") where 
h*^") is an integer-valued vector field. At this point the 
theory reads 
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/oo 
- oo 1 - , , 



°° 7=1 h(-') 



_J_^(A X h("))2 - ieA • I V A X h(") 

a—1 ' ' \a— 1 



f3 



(A X A) 



(B4) 



Again, we apply the Poisson summation formula Eq. ljB2|) and the integer fields h^") are replaced by continuous dual 
gauge fields h'"^ at the cost of introducing the term 27ri^^h'"^ • m'^"^ in the action and m^") are integer vortex 
fields. Then the partition function reads 



/OC ^ pOQ 

PAj] / ^ <5A.m(.),oexp(-5), 
-OO ^,— 1 — OO 

7—1 m(~0 



E ^ h^)^ - ieA • (^E A X h(") j + 2.1^: h(") • m(") + |(A x A)^ , 



(B5) 



where Sn,m is the Kronecker delta. The gauge symmetry of the integer valued fields h^"^ in Eq. ljB4|) must be 
preserved through this transformation. Hence, the action Eq. ljB5|) must be invariant under the gauge transformation 
h(") h*^"^ + Ax. The transformation produces the term J^r J2a SttIAx • m'"^ which must be zero to preserve gauge 
symmetry. A partial summation therefore produces the constraint A • m'"' = on each flavor of the vortex fields. At 
this point it is useful to write the theory in the Fourier representation and the action becomes 



(-TV N 

^ = E E ^si^^"^- X h^"^) • (Q- X h^l) + - E 



N 



le 
"2 



\a=l 



(B6) 



(Qq X Aq) • (Q_q X A_q) \ , 



where Qq is the Fourier representation of the lattice difference operator A. We choose the gauge A • A = and 

A • h*^"^ — which in the Fourier representation is Qq • Aq — and Qq • hq"-* = 0. We complete the squares in Aq 
and get the action 

( N N 

^ -E |E • • + - E [K' • ^'-t + • < 



N 



Aq-^ EQ.Xh("Ml?ql 



N 



(B7) 



N 



N 



T EQ.Xh(") I?qM5:Q_qXh 



(a) 

q 



where Dq = /3Qq • Q_q/2. After performing the gaussian integral in Aq, the action reads 

( N N 



(a) 

q 



N 



N 



\a=l 



\a=l 



At this point it is useful to introduce matrices and vectors in flavor indices. We write the action as 

•5 = E [^q Gq^f-q + iTTM^H^q + inH^M^q] , 



(B8) 



(B9) 
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where = {hq \ h.q \ . . . , hq^') and = (niq^'', niq'^'', . . . , niq"'') are vectors in flavor indices, and the matrix Gq 
is given by 



,(1) „,(2) 



2/3lii(i)|2 20 



Gq — 



2/3|V(i'P 2/3 

£_ 

2/3 2/3|V(2) 



2/3 



iQ<,r 



2/3 



V 



2/3 

el 

2/3 



2/3 
2/3 



2/3 
2/3 



2/3|V("-i>P 2/3 



2/3 
2/3 

el 

2/3 



iQq 



(BIO) 



2/3 2/3|V(")P 2/3 / 

where |QqP = Qq • Q-q- In flavor indices 77, A we write the matrix as Gq''"^' = S^'^Sri^x + c, where d*^''^ = „,|'^'(j,) 



and c = 1^ . We complete the square and obtain the expression 



[(^q + i7rMTGql)Gq(if_q + iTrA/^qGql) + TT^ Af ^ G" ^ M_q] 

q 

I 



(Bll) 



where is the solution of GqG"^ = / and / is the as ofF-diagonal elements (rj ^ A), where Y[a=0 = 

A'^ X iV unit matrix. Gaussian integration of i?q yields an 1. By dividing the numerator and denominator by 
action expressed in vortex variables of different flavors 

S^Y^ Tr^Af^Gq iM_q. (B12) 



Jl^^j^ d^"^ , we obtain the diagonal elements 



(GqO 



The matrix Gq ^ is found to have the following diagonal 
elements 



1 + -4 



(B15) 



(Gq 

and 



(Gql)('''^) 



N 



(B13) 



(B14) 



and the off-diagonal elements 



(G^l)(..A) ^ 



where rj ^ \. In total, the matrix (Gq ^)''''^'*'-' reads 



(B16) 



(G< 



Inserting the expressions for S-"^^ and c and multiplying by IQql"^ in the denominator and numerator, we obtain 



(B17) 



We introduce ~ E^=i ^-nd split the expression by partial fractioning and obtain the matrix 



1(2 



(B18) 



(B19) 



IQqP |QqP+e2vl/2_ 

This is the vortex interaction matrix given in Eq. \i'2,\ . Inserting this into Eq. (jB12|) . the partition function of the 
system is 



Z =W ^ (5A.m(°),oexp(-S'v), 



N 



q r),A=l 



IQqP 



|QqP + e2*2 



(B20) 
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which gives the action in the partition function in Eq. The above vortex action may be written in terms of 

charged and neutral vortex modes in a manner analogous to that of Eq. HA7|I . as follows 

27r^f3/^^ |QqP+m2 +^ 2|Qq|2 

Here, uiq = e^^*^. While the first, screened, term in general is present for all > 1, it is clear from the above 
formulation that the second, unscreened, term is only present provided N > 2. The factor 2 in the denominator in the 
unscreened terms is essential in order for the interaction terms between different vortex species to cancel out when 
ml = 0. 



APPENDIX C: GAUGE FIELD CORRELATOR 



Introducing Fourier-transformed variables in Eq. prior to integrating out the gauge field A, and adding source 

terms in the form of electric currents J coupling linearly to A, we obtain the action (note that for convenience, we 
have redefined the A-field in this Appendix by a factor as follows: A — > A = ^/PA, and then renaming A ^ A. 
Thus, here we also redefine the charge as follows: e — s- e = e/-\//3 which we then rename e e. In the end result we 
reinstate the original definitions. We thus have the action Eq. pOfl on the form 



= ^jinwh(") . hL") + TTi [m(") . hLt + • h(")l - f [a, . (Q_, X hLt) + (Q, X h(")) . A_ 



-\1 -A +1 .Al + l^^^A -A 



1 



(CI) 



Summations over indices (a, (3) G [1, . . . , N] is understood. Integrating out the gauge fields Aq, we get the action on 
the following form 



N 



^■^ - EE 



q a— 1 



2^|^(q)|2 q "-q^"^["M "-qT^"^-q 



Dq D q 

2|QqP 



(C2) 



where we have defined Dq = Jq + ieQq x hq and hq — X]a=i '^q ■ Thus, the last term in Eq. HC2(I may be written 

1 



Dq D q 



2|Qq 12"^"'' • Jq + hq • A_q + h_q • Aq - — hq • h_q. 



(C3) 



«q| ^|v<q 

where we have defined Aq = 2|^~p'Qq ^ "^q- Thus, the action may be written on the form 



2^|^(a)|2 q "-q^"^["^q "-q "^-q "q 



l^hq • A_q + h_q • AqJ + —hq ■ h_q 



2|Qq|2 



(C4) 



In Eq. ljC4p . a summation over indices (a,/9) is understood. We can now integrate out the dual gauge fields hq"'' to 
obtain 



where 



— Ila=l Em(°) '^A-M(°),0 exp [-5* Jeff] , 

Sj.^ = 5:[.^m(")5(-'')(q)mL^)-FA(j("),jLt) 



(C5) 
(C6) 



Here, we have introduced 



FA({Jq,J-q}) ^ \1 .7'^ - TTI 



2|Qq|2 



mL"^5("^'3) (q)L(f ) + lL^ (q)m(f )] + lL^ (q)L(f ) . (C7) 
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Here, the upper index denotes a "flavor" index (a,/?) € [1, . . . , A^] indicating which matter field species the fields 
above correspond to, and we have introduced the vector Lq"'' = Bq, a G [1, . . . , iV]. Using this particular property of 
Lq"-*, we may simplify FA({Jq, J_q}) somewhat, to obtain 



7M pA*^ 7!^ 

i^A({Jq,J-q}) ^ \1 -'' -TTl 



2|QqP 

In Eq. (|C8() . we have furthermore introduced 



mL"^ • Aq + A_q • m(")q F^") + A_q • Aq5. 



(C8) 



N 



/3=1 
N 

0,(1=1 



IQqP+mg' 



IQq 



(C9) 



The last equalities in Eqs. (|C9|I are found from using the definition of D'-'^'^\q) given in Eq. ((T^ . along with the 
definition of tp"^ given immediately after Eq. (|12|l . We have also introduced the transverse projection operator 



1/ QqQ-q 
IQqP ' 



(CIO) 



which appears due to the transversality of the currents J. Before doing fimctional derivations on FA({Jq, J-q}) it is 
useful to multiply out the term A_q • Aq and explicitly use the constraint V • J = before derivation. We find 



A_q • Aq 



^av>^ ^apr, njV r\P J>- 7') 

4|Qq|' 



4|Q, 14 



_ 7M pA"' jv 



Cql 4|Qq| 

In the cross terms between vortex fields and current fields there is no need to introduce the transverse projection 
operator, since the inner product automatically projects out the transverse part of J since the vortices form closed 
loops. 

We may now express the gauge field correlators formally 



J-a=Ja=0 



A-M(°).0 



a=l M(°) 



SFa SFa , 5^Fa 



exp [-^ooff] 



J-q=Jq=0 



where ^ocff is Eq. HC6|I with all currents set to zero. We get for the required functional derivatives 

"'iql 



^T^-^-q _ • m("V(") + 2A • ^5 



2|QqP ' • 



Multiplying out everything, after setting the currents to zero, we obtain 

- ^ [2 - ^'S] + Uv^^"V('^)«")777t^), 

where a summation over (a, /3) € [1, . . . , TV] is understood. Setting v = fi and summing over v, we get 

1 



IQqP 



[2 - e^S] 



IQqP 



(Cll) 



(C12) 



(C13) 



(C14) 



Using the definitions of U^"-* and S introduced above, and moreover reintroducing the original gauge field A and the 
original charge e (see start of the Appendix) , this becomes 

2 4/3^V |^(a)|2|^(/3)|2 



(Aq • A_q) 



/3 I IQqIM 

|QqP + m2 



■ IQqP (|QqP+m2)2 

2/3^2g2 G(+)(q) 



mi ' • m 



IQqP IQqP- 



(C15) 
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when we introduce 

GW(q) = m(")) . iV^'^f mL^.^)), (C16) 

(a) (3 

which are just Eqs. (|24|l and (|25|l . Note how this has the structure of a "Dyson's equation", where G^~^\q) plays the 
role of a matter field loop or " polarizability" , where the strength of the vertex is given by 2/37re/|Qq|. As we shall see 
below, a similar statement holds for the dual gauge field correlators, but there the vertex is a scalar of strength 7ri. 
The difference lies in the fact that while the dual gauge fields couple linearly to the vortex fields, it is the curl of the 
gauge field A that couples indirectly to the vortex fields (via the dual gauge fields) . 



APPENDIX D: DUAL GAUGE FIELD CORRELATORS 



The computation of correlators for the dual gauge fields h^") proceeds along the same lines as for the gauge field 
A, but sufficiently many details are different so we include it here for completeness. Introducing Fourier-transformed 
variables in Eq. I|10() . and after having added source terms in order to be able to compute correlators, the theory may 
be written on the following form 



2/3|^(a)|2 q -q 



ml") • hL"2 + mL"2 • hi") 



+ 2 



Dl) 



In Eq. ljDl|) . summation over indices (a, /?) G [1, . . . ,N] is understood and the last term appears after having integrated 
out the gauge field A. Notice how we in this case have added N source currents J^"), one for each vortex field m*^"). 
Integrating out the dual gauge fields hq") we get the action on the following form 



(D2) 



where we have defined 



(D3) 



Here, as in Appendix [HI the currents are divergence-free, V • J^") — 0,a € [1, . ■ . ,N], and the interaction matrix 
£)(",/3)(q) is defined in Eq. (O. As was the case for the A field correlator, the constraints on the currents J^") must 
be carefully kept track of when performing the necessary functional derivations in order to obtain the correlation 
functions. The generating functional is given by 



^■r — Ila=l Sm(°) '^A-M(°),0 exp [-S'jcff] 

Applying Eq. l(nTT|) to Eqs. and we find 

1 5^Zj 



(D4) 



1 ^ 

IT /J '5a-m(°),o 
^0 — 



a = l M(°) 



5Fh SFh 



exp [-S'ojoff] , 



(D5) 



-I J 



where SqjcS is Eq. ljD2p in the absence of source terms. We obtain from Eq. (jD3|l 

= i5(-«(q)J,-(^)(q)-.ri5(-«(qK(«, 



(D6) 



Here, we must keep track of two indices on the "magnetic" currents J^"), since we have one current (source term) 
coupling to each of the N dual gauge fields h^"), a € [1, . . . , A^]. The notation we use is that Jq^"-* is the ^ Cartesian 
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component of the current coupling to dual gauge field of flavor a[l, . . . , TV]. We use a corresponding notation for mq 
Moreover, PJ^'^ again is the transverse projection operator defined in Eq. (|C10|) . appearing due to the transversality 
of the currents J^"^. Multiplying all of this together, we find 

- ^D^^-^HcOPJ^" - 7r25("''')(q)i5('3''^)(q)(TO^('')m':'^^). (D7) 
Moreover, setting v — fi and summing over Cartesian components of the dual gauge fields, we find 

(h(j") • h^^l) = 5("''^)(q) - 7r25("''')(q)i5('3''^)(q)(mL''^ • in[^^), (D8) 

where we have used the fact that the trace of the projection operator is equal to 2. In Eqs. (|D7p and IjDSp , a summation 
over the indices (k, r/) € [1, . . . , A^] is understood. Note how this, as for the A-field correlator, has the structure of 
a "Dyson's equation", where the vortex correlator plays a role analogous to a matter field loop, or " polarizability" , 
with a scalar vertex (charge) of strength 7ri. This is simply a reflection of the fact that dual gauge fields h^"^ couple 
linearly to the vortex fields m^") ("magnetic currents"). The factor i in the strength of the vertex appearing in Eqs. 
(|10|l and ljDl|) gives rise to an "anti"-Biot Savart law between vortex segments mediated by the dual gauge fields. 



APPENDIX E: GENERALIZATION OF EQ. ^ IN THE PRESENCE OF INTER FLAVOR JOSEPHSON 

COUPLING 



In this Appendix, we consider the A^-flavor London superconductor model Eq. Q in the dual representation, 
including Josephson couplings between matter fields of different flavors. The Josephson coupling between B^"") and 
6'('') is local in space-time, represented in the Euclidean action by the terms g^°''^'> cos[6'("^ (r) — 9^"^^ (r)]. With N matter 
fields there will be N{N — l)/2 such terms. (Note how there are no such terms when N — 1). However, since these 
terms act as ferromagnetic couplings between the phase fields of different flavors, the critical properties of the model 
are preserved if we only include the terms that are "nearest neighbors" in flavor indices. This is precisely analogous to 
including only nearest- neighbor Josephson coupling in a Josephson junction array, but where the "lattice sites" now 
are represented by flavor indices. In this case, we have — 1 Josephson terms. Therefore, we consider the action 



N 



S^-Y^ /3| V'^"^ P cos(A6'(") - eA) 



a=l 



N-1 

E 



(v), 



3(6»('7)_6)(';+i)) + ^(Ax A)2, 



(El) 



where g^"^^ is the Josephson coupling. In the Villain approximation the model reads 



N-l 



fvAll f E n E ^M-s), 

^=E 



a=l n<") ')=! m<i 



J2 ^^^(A0(") - eA + 27rn("))2 + 'j^ ^(^'^^ - e^^+'^ + 2T,m^^^f + §(A x A) 



.a=l 



(E2) 



Here, n^") are integer vector flelds where a € [1, . . . , A^] and m^''^ are integer scalar fields with j] G [1, . . . , A^ — 1] 
which take care of 2n periodicity. We introduce the Hubbard Stratonovich fields v*^"^ and g'^''^ , and apply the Poisson 
summation formula so that they become integer fields 



^=E E 



a = \ v(a) i;=l q(^) 

N r . (a)^2 



^^+iv(").(A.(")-eA) 



r,= l 



(E3) 



^(A X A)^ 



At this point we organize the phase fields and perform This gives the following constraints on the integer fields 
partial summations so that they can be integrated out. 

A.v«=gW, 

A • v^**) q^**) - (E4) 
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To enforce these constraints we introduce the non- 
compact gauge fields h*^") and the integer fields B^^'^ 
where a € [1, . . . , A^] and rj £ [1, . . . , — 1] such that 



vW=BW+Axh«, 
v(';) ^ b(''' - b(''"1) + a X h(''\ 
vW = _B(^-i) + AxhW, 



(E5) 



where 77 e [2, . . . , — 1], and q^''^ = A • B^**' are instan- 
tons. These expressions may be simplified by introducing 



J 



the dummy fields = g(^) = and B^^) = B(^) = 0, 
so that the constraints become slightly more symmetric, 
given by 



A • v(") 



9 



B("-i) + A X h("\ 



(E6) 



where a € [1, . . . , N]. Including inter-fiavor Josephson 
couplings beyond "nearest- neighbor" would merely have 
led to redundant additional constraints in the problem. 
Expressed in the new fields, the partition function reads 



N-l 



/ n E n E ^-p(-^) 

a=l h(°) »?=1 Bl") 

i (Axh(°)+B(" 



B(")=B(" 

B(«-i))2 



N 



(3 



(E7) 



a = l 



(A X A)= 



The appearance of instantons and effectively compact dual gauge fields in the dual description when Josephson 
couplings are included, serves to illustrate what a non-perturbative effect this is. Instantons are singular objects and 
cannot possibly be introduced perturbatively. Moreover, once instantons are required, a compactification of the dual 
gauge fields is also required, a highly non-perturbative step. Therefore, when we consider models where Josephson 
coupling is absent, it is essential to be able to rule out completely inter-fiavor Josephson coupling on symmetry 
grounds a priori, and at the level of the bare action. (In systems with multicomponent electronic condensates, with 
N = 2 and where a weak inter-fiavor Josephson coupling must be included, it may be possible to see the resemblance 
of one phase transition for a neutral mode and one phase transition for a charged mode, such as we have presented 
for the zero- Josephson coupling case. This would be so in small enough systems with linear extent smaller than 
the Josephson length given by Aj"'' = which in this context may be viewed as setting the scale of 

the inter-instanton separation. This would be a finite-size effect. For bulk systems, the apparent neutral mode will 
eventually be suppressed, leaving one phase transition in the universality class of the inverted SDxy model.) 

We proceed by integrating out the original gauge field and apply the Poisson summation formula to introduce the 
integer vortex fields m^") and the integer fields where a £ [1, . . . , A^] and rj £ [1, . . . , A^ — 1]. The resulting theory 
is the generalization of Eq. (|10|l 



N N-1 



B(0)=B(«)=0 
2 N N-1 



. V ^ V (Axh(")+B(")^B("-i))^ " (A • BC'))^ 

r l,a=l / I T- I a=l r]=l 



(E8) 



(N N-1 
a — 1 77— 1 



iv) 



First we note that like in Eq. Hl()(l the integration of 
the gauge field A makes the algebraic sum of the dual 
gauge fields massive. This refiects the fact that the gauge 
field A cannot screen instantons, it can only screen vor- 
tices. Furthermore, in the limit g^"^") — > for all 77 £ 
[1, . . . , A^ — 1] there are no Josephson coupling terms and 
each field B^**^ is constrained locally so that A • B^''^ — 0. 
The representation B^**^ ^ A x b^**^ takes care of the con- 



straint, and the substitution h(") + b(") - b("-i) ^ h(") 
reduces Eq. (|E8p to Eq. (|10|l . Finally we consider the 
uncharged case, e — > for which it is useful to return to 
Eq. IIE3|) . integrate out the phase fields, and write the 
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theory in terms of the integer fields v^") and q^''^ 

N N-1 

^ = n E n E ^A-v(°).,(°)-9(— 1) exp(-5), 

a—l v(") '^7—1 gC^') 



TV 



E 



(E9) 



where is the Kronecker delta. We sum over the fields 
q'^''^ for all 77 S [1, . . . , — 1] and are left with the parti- 
tion function 



JV 



^ = n E A.v(°).o expl-S"), 



Q = l v(°) 

^=E E 



E 



(ElO) 



This is the theory of A'" current fields v*^") which individ- 
ually can form closed loops, or dumbbells starting and 
ending on instantons. There is only one remaining con- 
straint on the N matter fields in the problem, after the 
A'^ — 1 instantons have been summed out. (If we had had 
M < — 1 Josephson coupling between flavors to begin 
with, we would have had N — M remaining constraints 
in the problem after summing out the instantons). The 
one remaining constraint leaves only one phase transi- 
tion in the problem in the universality class of the SDxy 
model, in contrast to the A^ phase transitions we have in 
the complete absence of inter- flavor Josephson couplings. 
The local constraint , A-v^"' = forces each dumb- 
bell to form a closed loop with one or more dumbbells 
of any flavor. The instantons have been summed out of 
the problem, leaving as their only trace the possibility 
of having supercurrents change flavor, precisely what the 
inter-flavor Josephson coupling does when expressed in 
terms of the fields describing the supercurrents! 

It is the possibility of being able to " chop" closed cur- 
rent loops of a given flavor into dumbbell pieces start- 
ing and ending on instantons, and then joining together 
dumbbell configurations of different flavors, that facili- 
tates this. (It has already been noted^ that the dual field 
theory of a gauge theory with two complex scalar mat- 
ter fields minimally coupled to a compact gauge field, 
features two complex scalar dual matter fields coupled 
to one non-compact dual gauge field and an inter-flavor 
Josephson coupling between the two matter field compo- 
nents. This is the reverse of what we have shown in this 
Appendix for the case N — 2, and nicely demonstrates 
that "duality squared equals unity".) 

Taking the limit 5^''' ^ for all 77 S [1, . . . , A^- 1] con- 
straints v^") to be divergence free for all a e [1, . . . , A^], 
and the model thus reverts back to the loop-gas repre- 
sentation of A^ decoupled 3Dxy models. 



APPENDIX F: KOSTERLITZ-THOULESS 

TRANSITIONS IN N-FLAVOR 
SUPERCONDUCTOR IN TWO SPATIAL 
DIMENSIONS AT FINITE TEMPERATURES 



In 2 + 1 dimensions at finite temperature, the clas- 
sical critical behavior of the N-flavor superconductor is 
very different from the true 2 -I- 1-dimensional case, i.e. 
the quantum critical behavior taking place in two spa- 
tial dimensions at zero temperature. Let us first re- 
call some features of planar superconductivity. It is 
well known that in two dimensional models with a U{1) 
gauge symmetry there is no quasi-long range order at 
any finite temperatures because a gauge field coupling 
makes the interaction between topological defects expo- 
nentially screened^. The situation is however differ- 
ent if one takes into account the "out of plane" mag- 
netic field. That is, taking into account a third dimen- 
sion, a vortex in a thin superconducting film produces 
a "mushroom" -like magnetic field outside the plane, 
which as shown by Pearl^^ gives rise to a logarith- 
mic inter-vortex interaction at distances smaller than 
[penetration length]^ /[film thickness], while at a dis- 
tance larger than that the vortices interact via a l/r- 
law. Thus, for a thin film with a penetration length 
which is significantly larger than the sample size, a 
Kosterlitz-Thouless (KT) transition/crossover should be 
observableSSia. The same effect is also the reason for 
the appearance of vortices with long-range interactions 
in layered systems making them being essentially cou- 
pled U{1) models, where various KT transitions and 
crossovers were studied in numerous worksiSi^iiS. 

Here, we are interested in KT phase transitions in the 
N-flavor London superconductor in 2 + 1-dimensions in 
the regime when (i) the effect of "out of plane" Pearl field 
can be neglected (short penetration length limit or alter- 
natively a planar field theory without a third dimension), 
and when (ii) all components have different stiffnesses. In 
Ref. 1^ the case A^ = 2 in such a regime was considered. 
It was shown that the system has a KT transition into a 
state where quasi-long range order is established only in 
phase difference which produces a quasi-superfluid state. 
This state, however, is principally different from e.g. SSF 
to MSF transition considered in-^° and in this paper, 
because (i) the quasi-superfluid state considered in'^^ is 
a purely two-dimensional phenomenon (ii) in this state 
there is no true off-diagonal long-range order, (Hi) there 
is no phase transitions from superconductivity to super- 
fluidity in a planar system at finite temperature^. 

Let us now consider Eq. Q for the case A^ = 3, when 
I <c <C In the most interesting case 

of finite penetration length, the charged mode formally 
can never develop quasi-long range order. That is be- 
cause the composite single-quantum vortices (Afl^^^ — 
27r, A6'(2) = 27r,A6'(3) = 2n) and {A9^^'^ = -27r,A6l(2) ^ 
— 27r, AO'^^^ ~ — 27r) have finite energy and have only 
screened short range interaction. Thus, in the limit where 
the magnetic penetration length is short, such vortices 
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are always unbound at any finite temperature. We do not 
consider here the possibihty of a "would be" KT crossover 
which is possible in a charged system with significantly 
large penetration length^^. The absence of superconduc- 
tivity means that individually all phases are disordered 
and the system is not superconducting. However, con- 
sidering quasi-long range order in phase differences in 
d = 2, several interesting possibilities arise. Composite 
one flux quantum vortices have short range interactions. 
On the other hand, vortices with windings only in one or 
two phases excite neutral modes and thus can undergo a 
true KT transition. This opens up the possibility for a 
KT phase transitions associated with establishing quasi- 
long-range order in phase differences^i. The key feature 
of the N > 2 system where bare stiffnesses are different, 
is that, as discussed in Appendix IXI the neutral modes 
have also different stiffnesses [see Eqs. (15^ . (jS^ ]. 

Let us consider first the low temperature regime. 
Then the vortices with short-ranged interactions, namely 
(A6l(i) = 27r,A6l(2) = 27r,A6'(3) = 27r) and (A6l(i) = 
-27r,A6'(2) = -27r,A6'(3) = ~2n) are liberated, while 
vortices with phase windings only in one or two phases 
are bound into pairs of vortices and anti- vortices. In this 
state there is quasi-long range order in the phase differ- 
ences 6'^^^ ~ e^^\ 61(1) -e'ts), and 0(2) _ 5,(3)^ Recall that the 
gradient terms of neutral modes which follow from sep- 
arating of variable in GL functional and dropping terms 
describing charged modes in Eq. H53|l are 



1 |V'(i)p|V'(3)|2 

2' 



\I/2 

1 |V'(2)|2|V'(3)|2 

2 ^2 



(V(0(l)-0(3)))2 
~0(3)))2. 
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(Fl) 



From this, we have the following stiffnesses of neutral 
modes 



= 



|^,(1)|2|^(2)|2 
^2 ■ 

|^(1)|2|^(3)|2 
^2 ■ 

|^(2)|2|^(3)|2 



(F2) 



Consider now what will happen as the temperature is 
increased. Pictorially, we may illustrate this by con- 
sidering Fig. 1^ by slicing through the pictures at a 
given coordinate along the vortex lines, considering typ- 
ical cross sections. Upon increasing the temperature, 
first there will take place a deconfinement of vortex 
pairs (A6'(i) = 27r,A6'(2) 0,A6'(3) = 0)+ (Ae'(i) = 
— 27r, A^^(2) — 0, A0(^) — 0) because vortices in such a 
pair are bound by the two weakest neutral modes J^^ and 
[ji2 < ji3 < 1^(1) |2/2]. This wiU be accompanied 
by partial decomposition of deconfined thermally created 
composite vortices (A6'(i) = 27r, A6'(2) 2tt,A9^^^ = 
2tt) ^ (Ae'(i' = 27r, A6'(2) = 0, A6i(3) = 0) -I- (A6i(i) = 
0,A6'(2) = 27r,A6'(3) = 2tt), because a vortex {A9^^^ = 



0, A6l(2) = 27r, A6'(3) = 2tt) has the same neutral vor- 
ticity as a vortex (AO'-^^ = -2tt, A6'(2) = 0, A6'(3) = 0) 
namely — 27r windings in neutral modes 0(1-' — 0(2) and 
Qi^) _ 0(3) ^ This transition takes place at 



(1) _ ![rjl 
KT — 2'- 



2 



/,(2)|2 + |^(3)|2] (p3) 



This phase transition disorders the variable 0(i) and cor- 
respondingly eliminates quasi-long-range order in phase 
differences 9^^^ — 0(2) and 0(i) — 0('^). Consequently 
above Tj^^r^ only surviving neutral mode is associated 
with 0(2) — 0('^) . The remaining phase transition can be 
mapped onto that in = 2 system^i. Thus the second 
phase transition takes place at 



^KT - 2 ^2 



1^(3) I 



(F4) 



We can now solve the general problem of KT transitions 
in a system of N planar condensates with all different 
bare stiffnesses. In the general case of A^-flavor London 
model the temperatures of the lowest KT transition is 
given by 



(1) 



N 



KT 



^ |7/;(i)|2 
2 ^ 



N 



^2 



)|2 



a=2 



1)12 
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The subsequent KT transitions at higher temperatures 
are mapped onto the A^ — 1, iV — 2, . . . cases. Taking the 
A^ ^ oo hmit in Eq. ((F5|| and provided that |V'(i) P < \l/2 
we obtain 



(1)|2 



(F6) 



This expression quite remarkably shows that in the limit 
A^ — > oo, even in the system with short penetration 
length, Tj^x tends to the value in a neutral system with 
the bare stiffness IV'^^-'I- In contrast in the one compo- 
nent case with short penetration length the system does 
not exhibit a KT transition. 

In conclusion, we note that the KT transitions consid- 
ered in this Appendix are still significantly simpler than 
the situation arising in this model in three dimensions 
because, as we have considered in previous sections, in 
three dimensions the charged mode plays an extremely 
important role. It is precisely the interplay between neu- 
tral and charged modes which is particularly important 
in three dimensions and which gives the model a vari- 
ety of different phases and phase transitions. Also, we 
note that this situation is quite different from KT tran- 
sitions that are known to exist in the 2 -I- 1-dimensional 
N-component Chiral Gross-Neveu model^S where there 
is only one KT transition which occurs at finite tem- 
perature (when the system is effectively two-dimensional 
through dimensional compactification). 
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